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ABSTRACT 

The light curves of Type la supernovae (SNe la) are powered by the radioactive decay of 56 Ni 
to 56 Co at early times, and the decay of 56 Co to 56 Fe from ~ 60 days after explosion. We 
examine the evolution of the [Co III] A5893 emission complex during the nebular phase for 
SNe la with multiple nebular spectra and show that the line flux follows the square of the 
mass of 56 Co as a function of time. This result indicates both efficient local energy deposition 
from positrons produced in 56 Co decay, and long-term stability of the ionization state of the 
nebula. We compile 77 nebular spectra of 25 SN la from the literature and present 17 new 
nebular spectra of 7 SNe la, including SN 2014J. From these we measure the flux in the 
[Co III] A5893 line and remove its well-behaved time dependence to infer the initial mass 
of 56 Ni (M m ) produced in the explosion. We then examine 56 Ni yields for different SN la 
ejected masses ( M e j - calculated using the relation between light curve width and ejected 
mass) and find the 56 Ni masses of SNe la fall into two regimes: for narrow light curves (low 
stretch s ~ 0.7-0.9), Mn% is clustered near Mjyn ~ O.4M 0 and shows a shallow increase 
as M e j increases from ~1-1.4M 0 ; at high stretch, M e j clusters at the Chandrasekhar mass 
(1.4M 0 ) while Mjyfi spans a broad range from 0.6 — 1.2M 0 . This could constitute evidence 
for two distinct SN la explosion mechanisms. 

Key words: supernovae: general 
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1 INTRODUCTION 

Type la supernovae (SNe la) were instrumental to the d iscov ¬ 
ery of the accelerating ex pansion of the Universe dRiess et al] 
1 19981 : IPerlmutter et al.1 19991) and remain key tools for character¬ 
izing the precise co s mology of the Universe dKessler et alJ 1 20091 : 
ISullivan et alj 1201 ll : iRest et ail 120141 : iBetoule et al.1 l2014ff Their 
cosmological utility is facilitated both by their intrinsic brightness 
{Mb ~ —19 at peak) and the relative uniformity of their peak 
brightnesses. More importantly, their luminosity div ersity i s tightly 
correlated with the width of the optical light curve dPhillipsIf1993h . 
The physical origin of this width-luminosity relation (WLR) has 
long been a subject of debate and is intimately tied to the progeni¬ 
tor system of SNe la and the physical mechanism that triggers the 
explosion. 

SNe la are widely believed to result from the thermonuclear 
disruption of a carbon-oxygen (CO) white dwarf dHovle & Fowled 
I960), which has recentl y been supported observationally for th e 
very nearby SN 201 lfe dBloom et all 120121: iNugent et all 1201 ll) . 
The CO-rich material in a white dwarf (WD) is supported against 
gravitational collapse by electron degeneracy pressure. A stable 
isolated WD lacks the internal pressure and temperature neces¬ 
sary to fuse CO to heavier elements (but see IChiosi et al.ll2015l) . 
In SNe la, this balance is upset by interaction with some binary 
companion, which triggers runaway nuclear fusion of the CO ma¬ 
terial to heavier elements, particularly iron group elements (IGEs) 
dominated by radioactive 56 Ni. The energy from fusion unbinds the 
star and ejects material at ~ 10 4 km s _1 . As the ejecta expand the 
decay of 56 Ni to 56 Co (with half-life of tj /2 = 6.08 days) releases 
energy into the ejecta which powers the optical lightcurve ofthe SN 
for the first few weeks after explosion ( [Colgate & Me Keel 196911 . in¬ 
cluding the luminous peak. At later epochs (t > 60 days past peak 
brightness), the SN la lightcurve is powered by 56 Co decay to 5(, Fe 
(with half-life of ti /2 = 77.2 days). Thus understanding the origin 
of the trigger mechanism and the amount of 56 Ni produced in the 
explosion would reveal the critical elements that make SNe la such 
excellent cosmological tools. 

The nature of the CO-WD binary companion is directly re¬ 
sponsible for the event that triggers the SN la exptosion. On e pos¬ 
sible scenario is the single degenerate (SD; Whelan & Ibeij l 19731 : 
iNomotdl 19821) scenario in which a CO-WD steadily accretes from 
a non-degenerate (main sequence or giant-like) companion until 
the central density of_the WD exceeds the critical density for car¬ 
bon ignition (e.g. jGasques et aPl2005l ) as the mass approaches the 
Chandrasekhar mass {Mw d 1 AMq ). In this scenario, the WLR 
has been proposed to arise from stochastic variations in the time 
at which the nuclear burning front within the exploding WD tran¬ 
sitions from sub-sonic to super-sonic^j- the so-called deflagration 
to detonation transi t ion (DDT; e. g., B1 inni kov & Khokhlov JV8(i 
Ropke & Niemevetl 20071 : IKasen & Wooslevl 120071 : iKasen et al.1 


(Guillochon et al. 201C 

; Loren-Aguilar et al. 2010; Pakmor et al. 

20 ljl 20 ll; Moll et al. 

2014; Raskin et al. 2014). This scenario is 


20091 : Sim et al Jl20 1 3l) . Variations in the time of the DDT result in 


different amounts of 56 Ni being produced, yielding diff erent peak 
magnitudes and light curve widths for SNe la (though fsim et al.1 
I2013I . do not recover the observed WLR). 

The other popular scenario f or SN la progenitor sys- 
tems is the double degenerate (DD; Tutukov & Iungelson Jl97g; 

nkll984) 


iTutukov & Yungelsonll979l : llben & Tutukov! 1984IWebbink8 198. 

scenario in which two WDs in a close binary merge after or¬ 
bital decay due to gravitational radiation. Some recent simu¬ 
lation results have shown that a violent merger of the two 
WDs produces_“hot spots” whicfyexceed the critical temperature 
and density (Seitenzahl et al. 2009a) needed to ignite CO fusion 


inherently not tied to Mch, but instead could produce explosions 
with varying luminosities and light curve widths simply due to the 
variation in mass of the progenitor system iRuiter et al.i2013l) . Gen¬ 
erally for the DD scenario, the WD undergoes a complete detona¬ 
tion and th e amount of 56 N i produced depen ds on the mass of the 
progenitor dFink et all20ld:lsim et all2010l) . 

Finally, it is important to also consider the double detonation 
(DDet) mechanism for triggering the WD explosion. In this sce¬ 
nario, helium-rich material accreted onto the surface of the white 
dwarf (either from a He-rich main sequence or giant star or He- 
WD) could ignite and send a shockwave into the core of the star. 


core which initiates the thermonuclear runawai 

process dLivne 

1990 

Iben & Tutukov 1991; Wooslev & Weaver 

1994 Fink et al. 

2010 

Wooslev & Kasen 2011). This mechanism could arise from 


SD or DD systems, and is not tied to Mch- Additionally, this 


of high-velocity features in early SN la s 

Dectra iMazzali et al. 

2005 

; Maguire et al. 

203; 

Childress et al. 

2013d; Marion et al. 

2013 

; Childress et al 

20144 

Maguire et al. 2014; Pan et al. 2015a; 

Silverman et al. 2015 

). 


While much of the debate about SN la progenitors in 
the previous decade revolved around which single scenario 
was responsible for SNe la, recent results have pointed to¬ 
ward multiple progenitor channels being realized in nature. 
SN la rates studi es yiel ded evidence for both short - and long 


lived progenitors (Mannucci et al. 2005; Scannapieco & Bildsten 

2005: ISullivan et alj 120061 : iMannucci et al J 120061 : lAubourg et al 


2008). The lack of a detected companion star to the progeni- 
jLietal 201 i. and in SN la remna n ts (e g . 
'al 120121: iKerzendorf et al.l 120 1 1 120 H 12014T) 


tor of SN 201 lfe dLi et al 

[Schaefer & Pagnottal 120121 ; 
present individual cases where the DD scenario seems necessary, 
while stro ng emission from ci r cum-stellar material in some nearby 
SNe la dHamu v et al] 1 20031: lAldering et all 120061 : iDildav et~akl 
l2012l : lsilverman et al J2013c fbh seems to indicate clear cases of the 
SD scenario. 

For peculiar white dwarf supernovae, like the Type-lax SN 
2012Z, a luminous proge nitor syste m has been detected and inter¬ 
preted as the donor star dMcCullv et alj|2014al) . Similarly, shock 
interaction of SN ejecta with a (non-degenerate) companion star 
has been detected in the e arly light cu rve of another peculiar, low- 
velocity white dwarf SN dCao et alJl20o) . However such shock 
interaction is distinctly absent for several other SNe la observed 
continuo usly thr ough th e epoch of first light with the Kepler satel¬ 
lite dolling et alj|2015h . Additionally, a general dic hotomy in the 
spect roscopic properties of SNe la appears evident dMaguire et al.1 
120141) . Thus numerous lines of evidence now point to multiple 
SN la progenitor channels being active. 

Variations in progenitor masses between different explosion 
mechanisms will manifest as d i versity in the bolometric light 
curves of SNe la dAmettl 1 19821 : Ijeffervl [l99^: IStrit zinger et all 
l2006al : iRopke et alJl2012m Recently. IScalzo et al.l d2014al) demon- 
strated that the ejected mass - hence the progenitor mass - of 
a SN la could be recovered to 10-15% precision, as tested on 
bolometric light curves derived from radiative transfer modelling 
of SN la explosion models with known input progenitor_mass. 
Applyin g the same modelling technique to real data, IScalzo et al.l 
d2014al) found evidence that the ejected mass varies in the range 
0.9-1. 4 Mq among spectroscopically normal dBranch et al.ll 1993h 
SNe la and that the ejected mass also correlates strongly with the 
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SN la 56 Ni masses from nebular 56 Co emission 3 


light curve width parameter used to standardize SN la distances in 
cosmology. The correlation between ejected mass and light curve 
width was exploited by lScalzo et alj d2014bh to measure the SN la 
ejected mass distribution: they found that 25-50% of all normal 
SNe la eject sub-Chandrasekhar masses, with most of the rest 
being consistent with Chandrasekhar-mass events (this is consis¬ 
tent with constraints from Galacti c chemical evolution b ased on 
Mn/Fe in the solar neighborhood - ISeitenzahl et all2013ah . Super- 
Chandrasekhar-mass SNe la were found to be very rare, at most a 
few percent of all SNe la, consistent with previous measurements 
of the relative rate dScalzo et all201 j ). 

The diversity in ejected mass suggests a corresponding di¬ 
versity in explosion mechanisms among normal SNe la. Fur¬ 
ther information about the explosion mechanism may also_be en¬ 
coded in_the_ge^_absolute magnitude distribution I Ruit er et al.l 


120131: PiroetaL 201 4) . the diversity in early SN la light curves 
(DessartetaL 2014c ), or in the relation between 56 Ni and ejected 


mass 
56 


Isim et al . 2oT( A I Ruiter et alho i ^lScalzo et all2014al ). The 

Ni mass is most commonly inferred from the peak absolute mag¬ 
nitude of the supernova dAmettlll 982|) , although with some model- 
depen den t syste matic errors (Branc hi 1 19921: iHoeflich & Khokhlovl 
1 19961: iHowell et al.l l2009h . The 56 Ni mass can also be inferred 
from detailed modelling of photospheric phase spectral times se- 
ries 1 Stehleetal^ 2005Mazzalietal 120081 : iTanaka et al] [ 20 TTI; 
ISasdelli et al.l2014 ; Blondin et al.l2015 ). Reliable alternative meth¬ 
ods for measuring 5G Ni masses, with different model-dependent 
systematics, can thus in principle help to shed light on the explosion 
mechanisms and progenitor properties of SNe la. 

In this work, we show that the amount of 56 Ni produced in the 
SN la explosion can be measured directly from signatures of its de¬ 
cay product 56 Co in nebular phase spectra of SNe la. Specifically, 
we employ the flux of the [Co III] A5893 line in spectra of SNe la 
in the nebular phase (t > 150 days past maximum brightn ess) as 
a diagnostic of the mass of 56 Co at a given epoch. IKuchner et al .1 
d 19941) showed that the ratio of the flux of this line to the Fe III 
line at A4700 A as a function of SN phase followed the expected 
temporal evolution of the Co/Fe mass ratio, which they used as 
evidence for the presence of 56 Ni generated in the SN explosion. 
More recently the presence of S(, Ni has been directly confirmed 
through 7 -ray line emi ssion from 56 Ni dDiehl et alj 20141) and 56 Co 
dChurazov et al.ll2ol4l) lines observed by the INTEGRAL satellite 
for the very nearby SN 2014J. 

Previous studies of SN la nebular spectra have collected a 
modest sample of spectra (a few dozen) from whic h important sci¬ 
entific results were derived. iMazzali et alj dl998l) found a strong 
correlation between the width of nebular emission lines (specif¬ 
ically the Fe III 4700 feature) with the SN light curve stretch, 
constituting evidence for greater 56 Ni production in more lumi¬ 
nous slow-declining SNe la. This result was combined with de¬ 
tailed modelling of nebular spectra (especially the 7380 A nebular 
line presumed to arise from stable 58 Ni) to infer a common ex¬ 
plosion mechanism for SNe la jMazzali et al.l2007h . Nebular spec¬ 
tra have also been employed t o place upper limits on hydrogen in 
the vicinity of n ormal SNe la llL eonard 20071. IShappee et al.OOl 3l : 


Isllverman et alj|2013al : [Lundqvislet al 


201 a) - The lack of hydro¬ 


gen in normal SNe la is in contrast to the strong hydrogen lines 
found in late phase spectra of SNe la which exhibited_ strong in¬ 
teraction during the photospheric phase (Silverman e t al.ll2013bl ). 
Velocity shifts in the purported Ni 7380 A nebul ar line were 
used to infer asymmetry in the inner core of SNe la dMaeda et al.1 
l 2010 al lbh. which was also found to correlate with the optical colour 
and Si 6355 A velocity gradient during the photospheric phase 


dMaeda et alj|201 lh . These line velocity shifts were found to also 
corre late with photospheric phase spectropolarimetry (Mau nd et al.l 
1M3). indicating a general correlated asymmetric geometry for 
SNe la. These early results have generally been supported with 
greater statistics afforded by new large data sets such as the Cf A 
sample dBlondin et al]2012n and BSNIP dsilverman et alj2013ah . 

Until recently, the nebular line at 5890 A was not frequently 
emphasized as a diagnostic of 56 Co due to its presumed association 
with emission from sodium (Kuchner et al. 1994; McClelland et al. 


l 20 ll are noteworthy exceptions). However 


Dessart et al. 


1 2014bl) 


showed definitively that this line arises primarily from cobalt for 
the majority of SNe la. We exploit this result to use the [Co III] 
A5893 line as a diagnostic of Sb Ni from a large sample of nebular 
SN la spectra compiled from both new observations and from the 
literature. Equipped with a sample of 77 spectra of 25 SNe la from 
the literature and 17 new spectra of 7 SNe la, we calculate the ab¬ 
solute flux of the nebular [Co III] A5893 line by scaling the spectra 
to flux-calibrated photometry measurements. With these calibrated 
fluxes we show that the temporal evolution of the absolute [Co III] 
A5893 line flux is highly consistent for SNe la with multiple nebu¬ 
lar spectra. We exploit this result to place measurements from dis¬ 
parate epochs on a common scale. This allows us to meaningfully 
compare the line fluxes in order to determine the relative amount of 
S6 Ni produced by each SN la in our sample. 

In Section [2] we present our compilation of literature SN la 
spectra and the new nebular SN la data released here. Section [3] 
presents our method for measuring the [Co III] A5893 flux from 
the spectra and scaling the spectra with the SN la photometry. We 
examine the temporal evolution of the [Co III] A5893 for SNe la 
with numerous nebular observations in Section [4] We then infer 
56 Ni masses for our SN la sample in Section [5] and discuss the 
implications and limitations of our results in Section[ 6 ] Finally we 
conclude in Section[7] 


2 SN la NEBULAR SPECTROSCOPY DATA 

The analysis in this work relies on a compilation of SN la nebular 
spectra from the literature as well as new observations. The full 
sample of literature and new late phase spectra are presented in 
Table [I] 
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Table 1. New and Literature Late Phase SN la Spectra 


SN 

Phase t a 
(days) 

Obs. Date b 

Spec. Ref. c 

SN 1990N 

160 

19901217 

BSNIP 


186 

19910112 

Gomez & Lopez (1998) 


227 

19910222 

Gomez & Lopez (1998) 


255 

19910322 

Gomez & Lopez (1998) 


280 

19910416 

Gomez & Lopez (1998) 


333 

19910608 

Gomez & Lopez (1998) 

SN 1991T 

113 

19910819 

BSNIP 


186 

19911031 

BSNIP 


258 

19920111 

Gomez & Lopez (1998) 


320 

19920313 

BSNIP 


349 

19920411 

BSNIP 

SN 1994ae 

144 

19950422 

BSNIP 


153 

19950501 

CfA 

SN 1995D 

277 

19951124 

CfA 


285 

19951202 

CfA 

SN 1998aq 

211 

19981124 

Branch et al. (2003) 


231 

19981214 

Branch et al. (2003) 


241 

19981224 

Branch et al. (2003) 

SN 1998bu 

179 

19981114 

CfA 


190 

19981125 

CfA 


208 

19981213 

CfA 


217 

19981222 

CfA 


236 

19990110 

BSNIP 


243 

19990117 

CfA 


280 

19990223 

BSNIP 


329 

19990413 

Cappellaro et al. (2001) 


340 

19990424 

BSNIP 

SN 1999aa 

256 

19991109 

BSNIP 


282 

19991205 

BSNIP 

SN 2002cs 

174 

20021106 

BSNIP 

SN 2002dj 

222 

20030201 

Pignata et al. (2008) 


275 

20030326 

Pienata et al. (2008) 

SN 2002er 

216 

20030410 

Kotak et al. (2005) 

SN 2002fk 

150 

20030227 

BSNIP 

SN 2003du 

109 

20030823 

Stanishev et al. (2007) 


138 

20030921 

Anupama et al. (2005) 


139 

20030922 

Anupama et al. (2005) 


142 

20030925 

Stanishev et al. (2007) 


209 

20031201 

Stanishev et al. (2007) 


221 

20031213 

Stanishev et al. (2007) 


272 

20040202 

Stanishev et al. (2007) 


377 

20040517 

Stanishev et al. (2007) 

SN 2003hv 

113 

20031228 

Leloudas et al. (2009) 


145 

20040129 

Leloudas et al. (2009) 


323 

20040725 

Leloudas et al. (2009) 

SN 2004bv 

171 

20041114 

BSNIP 

SN 2004eo 

228 

20050516 

Pastorello et al. (2007) 

SN 2005cf 

319 

20060427 

Wang et al. (2009) 

SN 2007af 

103 

20070620 

CfA 


108 

20070625 

CfA 


120 

20070707 

BSNIP 


123 

20070710 

CfA 


128 

20070715 

BSNIP 


131 

20070718 

CfA 


151 

20070807 

BSNIP 


165 

20070821 

BSNIP 


308 

20080111 

CfA 
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Table 1 (cont’d) 


SN 

Phase t a 

Obs. Date b 

Spec. Ref. c 


(days) 



SN 2007gi 

161 

20080115 

Zhang et al. (2010) 

SN 20071e 

317 

20080827 

BSNIP 

SN 2007sr 

177 

20080623 

CfA 

SN 20091e 

324 

20101016 

T15b 

SN 201 lby 

206 

20111202 

Silverman et al. (2013a) 


310 

20120315 

Silverman et al. (2013a) 

SN 201 lfe 

74 

20111123 

Shappee et al. (2013) 


114 

20120102 

Shappee et al. (2013) 


196 

20120324 

Shappee et al. (2013) 


230 

20120427 

Shappee et al. (2013) 


276 

20120612 

Shappee et al. (2013) 


314 

20120720 

Taubenberger et al. (2015) 

SN 201 liv 

318 

20121024 

T15b 

SN 2012cg 

330 

20130507 

M15 


342 

20130513 

T15b 

SN 2012fr 

101 

20130221 

This work 


116 

20130308 

This work 


125 

20130317 

This work 


151 

20130412 

This work 


222 

20130622 

This work 


261 

20130731 

This work 


340 

20131018 

This work 


357 

20131103 

M15 


367 

20131114 

This work 

SN 2012hr 

283 

20131006 

This work 


368 

20131230 

This work 

SN 20 Baa 

137 

20130710 

This work 


185 

20130827 

This work 


202 

20130913 

This work 


342 

20140131 

This work 


358 

20140216 

M15 


430 

20140422 

M15 

SN 2013cs 

320 

20140322 

This work 

SN 2013dy 

333 

20140626 

Pan et al. (2015a) 


419 

20140920 

This work 

SN 2013gy 

276 

20140920 

This work 

SN 2014J 

231 

20140920 

This work 

Note. — “ 

With respect to date of 5-band peak brightness. 

° Observation dates that are 

italicized are 

not used to measure Mjvi, and 

are only employed in Section 

11 


c BSNIP: 

Silverman et al. 

(20123); 

CfA: Matheson et al. 1 20081); 

Blondin et al. 

(20li); M15 

Maguire et al. (2015), in preparation; 


T15b: Taubenberger et al. (2015b), in preparation. 


tion change (we examine potential impact from the latter effect in 
Section [6.2t . 

To begin compiling a sample that meets these phase cri teria, 
we performed a large query of the WISeRElf] dYaron & Gal-Yaml 
12013 ) database to search for SNe la with two spectroscopic obser¬ 
vations separated by at least 100 days - assuming the earlier one 
would be near maximum light, this singles out SNe la with nebu¬ 
lar spectra. We then require SNe to have photospheric-phase opti¬ 
cal light curves sufficient to robustly establish light curve stretch, 
colour, and the date of maximum light using SiFTO JConlev et al] 
I2OO8I ). We also require the spectra to have sufficiently high signal- 
to-noise so that the [Co III] A5893 line can be well fit using a Gaus¬ 
sian fitting procedure (see Section[3]i. SN 2006X was excluded (de¬ 
spite having numerous neb ular spectra ) due to significant variabil¬ 
ity in its sodium features j Patat et al. 2007 ) and a rather s ignifi- 
cant light echo JWang et~aD [2008al ; ICrotts & Yourdonll2008l) . both 
of which might affect the time evolution of the [Co III] A5893 flux. 

Finally, we excise any SNe la which are spectroscopically 
peculiar in the n ebular p hase: SNe la similar t o SN 1991bg 
dFilippenko et alJ 1 1992bl : iLeibundeut et al.1 1 19931) exhibit ex¬ 
tremely narrow Fe lines and unusual line ratios; Ia-CSM SNe 
dSilverman et al.l [2013d) are excluded due to possible impa ct of 
CSM on the nebular emission; SNe lax dFolev et al] 1201 31 ) are 
excised as these probably arise from a different physical mecha¬ 
nism than norm al S Ne la; candidate super-Chandrasekhar SNe la 
dHowell et al.l 120061) are excised due to their unusual nebular 
spectra 1 TaubenbergMetaL 20131). SN e la similar to SN 1991T 
dPhillips et al.l 1 19921 : iFilippenko et al.1 Il992al) or SN 1999aa are 
however included in the sample, as their ionization structure ap¬ 
pears to be similar to “normal” SNe la. 

In summary, the selection criteria for our sample of literature 
nebular SN la spectra are: 


• Phase (with respect to B-band maximum light) in the range 
+ 150 < t < +400 

• Well-sampled multi-colour photospheric phase light curve 
(such that the light curve fitter SiFTO converges) 

• Sufficient spectrum S/N to measure the [Co III] A5893 line 
center and width 

• No spectroscopic peculiarity, except SN 1991T-like 

The full sample of spectra which meet these criteria are presented 
in TableQ] and comprise 77 spectra of 25 SNe la from the literature. 

Finally we note that two of the SNe in our sa mple had promi- 
nent light ec hoes at late times; SN 19 91T dSchmidt et al.l 1994h and 
SN 1998bu dSpvromilio et alj|2004l) . For both of these SNe, the 
light echo contributions are negligible at the spectroscopic epochs 
we employ. 


2.1 Compilation of Literature Data 

For reasons outlined in Section [4] the earliest epochs from which 
we can use [Co III] A5893 line fluxes is at phase t = +150 days. 
In practice, we found for most spectra beyond t ~ +400 days that 
the [Co III] A5893 flu x was too weak to be usable for our preferred 
analysis. Furthermore. [Taubenberger et al.1 J20l^ l recently showed 
that the t = +1000 day spectrum of SN 201 lfe showed dramatic 
changes in its structure, likely arising from a change in the ioniza¬ 
tion condition of the nebula. Indeed, this ionization change appears 
evident in the t = +590 day spectrum presented in ?, and we see 
evidence for the onset of this change shortly after t « +400 days 
in the data gathered for this analysis. Thus we excise data later than 
t ~ +400 days as unreliable due to low signal and likely ioniza- 


2.2 New SN la Nebular Spectroscopy 

We obtained new late phase (+50 < t < +150 days) and neb¬ 
ular (t > +150 days) spectra of 7 nearby SNe la from numerous 
telescopes. These spectra have been released publicly on WISeREP, 
with several spectra of SN 2012fr released through PESSTO's ESO 
data releases Q Information about observation details are presented 
in Table [2] and a plot of the spectra is shown in Figure [T| We note 
these spectra have not been rescaled to match observed photometry. 


1 http://wiserep.weizmann.ac.il 

2 www.pessto.org 
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Wavelength (A) Wavelength (A) 


Figure 1. New late phase and nebular spectra of SNe la presented in this work. All spectra are publicly available on WISeREP (except the +356 day spectrum 
of SN 2012fr from K15). Some spectra have been slightly binned (to ~ 5 A) for visual clarity. 


Several late phase spectra of very nearby SNe la we r e col - 
lected with the Wide Field Spectrograph (WiFeS; iDopita et al.1 
I 200 II 2010 I) on the ANU 2.3m telescope at Siding Spring Observa¬ 
tory in northern New South Wales, Australia. Observations were 
performed with the B3000 and R3000 gratings with the RT560 
dichroic, giving wavelength range of 3500 A-9800 A, with reso¬ 
lution of 0.8 A and 1.2 A on the blue and red ar ms, r espectively. 
Data were reduced using the PyWiFeS package dChildress et al .1 
2014bl). and spectra we re extracted using our custom GUI (see e.g., 
Child ress et all [20 1 3cJ) . We generally observed during very dark 
nights (moon illumination less than 20%) when the seeing was fa¬ 
vorable (l'.'5-2'.'0). We note that the WiFeS spectra of SN 2012hr 
and SN 2013cs have too low signal-to-noise to obtain a reliable 
measurement of the [Co III] A5893 line flux, but we release them 
publicly (on WISeREP) here. 

New nebular spectra for three nearby SNe la were collected 
with DEIMOS d Faber et alj 120031) on the Keck-II telescope on 
Mauna Kea, Hawaii. Observations were conducted with a l'/5 
longslit, the 600 1/mm grating with a central wavelength of 5000 A 
and with the GG410 order blocking filter, yielding a wavelength 
range of 4000 A-7650 A with 0.6 A resolution. Datajwere re¬ 
duced using standard techniques in IRAF (see e.g. Jchildress et all 
l2013al ). with the blue and red chips reduced separately then com- 
bin ed as a final step . We employed the Mauna Kea extinction curve 
of iButon et alj d2013l) . Our observations come from a single night 


on Keck (2014-Sep-20 UTC) when conditions were less favor¬ 
able (high humidity and thick clouds, ultimately 50% time lost to 
weather) but with a median seeing of 0"9. 

Five additional late phase spectra of SN 2012fr were collected 
as part of the Public ESO Spectroscopic Survey of Transient Ob¬ 
jects (PESSTO: ISmartt et al.lbo 1 si) durin g early 2013, and reduced 
with the PESSTO pipeline as described in lsmartt et al J d2015h . One 
spectrum of SN 2012fr and two spectra of SN 2013aa were ob¬ 
tained in 2013 using the Robert Stobie Spectrograph on the South 
African Large Telescope (SALT), and reduced using a custom 
pipeline that incorporates PyRAF and PySALT dCrawford et al J 
l2010h . One spectru m of SN 2012hr was obtained with Gemini 
GMOS dHook et afll2004l) using the 0.75" longslit with the B600 
and R400 gratings in sequence to yield a spectral coverage from 
4000 - 9600 A, under program GS-2013B-Q-48 (PI: Graham) - 
the spectrum was reduced using the Gemini IRAF package. 


In the analysis below we also include nebular spectroscopy 
samples from forthcoming analyses by Maguire et al. (2015, in 
prep. - hereafter M15) and Taubenberger et al. (2015b, in prep. 
- hereafter T15b). The M15 sample were obtained jrver a multi¬ 
period program at the VLT using XShooter 1 VernetrtjdJ^OLl ), and 
were reduced with the XShooter pipeline ( Modigliani et al. ] F20ld) 
using standard procedures (as in lMaguire et al. 20131) . The T15b 
sample were observed as part of a separate multi-period program 
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Table 2. Observation details for new late phase SN la spectra 


SN 

Phase 

(days) 

Obs. 

Date 

Telescope 
/ Instrument 

SN 2012fr 

+51 

2013-Jan-02 

NTT-3.6m / EFOSC 


+62 

2013-Jan-13 

NTT-3.6m / EFOSC 


+77 

2013-Jan-28 

NTT-3.6m/ EFOSC 


+ 101 

2013-Feb-21 

NTT-3.6m/EFOSC 


+ 116 

2013-Mar-08 

ANU-2.3m / WiFeS 


+ 125 

2013-Mar-17 

NTT-3.6m/EFOSC 


+151 

2013-Apr-12 

ANU-2.3m / WiFeS 


+222 

2013-Jun-22 

ANU-2.3m / WiFeS 


+261 

2013-Jul-31 

ANU-2.3m / WiFeS 


+340 

2013-Oct-18 

SALT / RSS 


+367 

2013-Nov-14 

ANU-2.3m / WiFeS 

SN 2012hr 

+283 

2013-0ct-06 

Gemini / GMOS 


+368 

2013-Dec-30 

ANU-2.3m / WiFeS 

SN 20 Baa 

+ 137 

2013-Jul-10 

SALT / RSS 


+ 185 

2013-Aug-27 

SALT / RSS 


+202 

2013-Sep-13 

ANU-2.3m / WiFeS 


+342 

2014-Jan-31 

ANU-2.3m / WiFeS 

SN 2013cs 

+320 

2014-Mar-22 

ANU-2.3m / WiFeS 

SN 2013dy 

+419 

2014-Sep-20 

Keck-II / DEIMOS 

SN 2013gy 

+276 

2014-Sep-20 

Keck-II / DEIMOS 

SN 2014J 

+231 

2014-Sep-20 

Keck-II / DEIMOS 


using F0RS2 on the VLT, and data w ere reduced with standard p ro- 
cedures similar to those employed in lTaubenbereer et aI1 li2013l) . 


3 NEBULAR LINE FLUX MEASUREMENTS 

3.1 The [Co ill] A5893 line in the nebular phase: a radiative 
transfer perspective 

The current study was motivated by the disappearance of Co III 
lines in nebular time series spectra, most notably the feature near 
5900 A. Previous literature analyses have attributed this feature al¬ 
ternately to Co III and Na I, so we turned to radiative transfer cal¬ 
culations to settle this ambiguity. 

We employed the time-dep endent radiative transfer code 
CMFGEN jHillier & DessartlfeoiTl) . which solves the time depen¬ 
dent radiative transfer equation simultaneously with the kinetic 
equations. Given an initial explosion model, we self-consistently 
solve for the temperature structure, the ionization structure, and the 
non-LTE populations, beginning the calculations at 0.5 to 1 day af¬ 
ter the explosion. The modelling assumes homologous expansion, 
typically uses a 10% time step, and no changes are made to the 
ejecta structure (other than that required by homologous expansion) 
as the ejecta evolve in time. Further details about the general model 
set up, and model atoms, can be found in iDessart et alj J20l4d ). 
We deployed CMFGEN on a delayed-deto nation model (DDC10 - 
iBlondin et al.ll2013l : IDessart et al.ll2014ch at very late phases and 
examine the contribution of various ions to the nebular emission 
spectrum. Radiative transfer calculations for this model, and simi¬ 
lar models but with a different initial 56 Ni mass, have shown fa¬ 
vorable agreement with observations dBlondin et alj I20T3L120151 : 
IPessart etalj|2014(l blfal). 

In Figure[2]we show PPC10 modeled with CMFGEN at phases 
+ 126 days (left panels) and +300 days (right panels). The top pan¬ 
els in each column show the integrated PPC10 model flux com¬ 
pared to observations of nebular phase SNe la at similar phases, 


while the bottom panels show the line emission from individual 
ions (note this can exceed the integrated flux due to the net opac¬ 
ity encountered by photons following their initial emission). The 
+ 126 day model shows particularly good agreement with the data. 
At +300 days the model shows some discrepancy with the data, 
particularly in the ionization state of the nebula. 

Most importantly, the radiative transfer calculations show that 
the emission feature near 5900 A is clearly dominated by Co III 
emission, with little or no contamination from other species. Few 
other features in the optical region of the spectrum show such clean 
association with a single ion. 

For later aspects of our analysis we require the velocity cen¬ 
ter of the nebula, which we calculate from the [Co III] A5893 line. 
To do so requires an accurate calculation of the mean rest wave¬ 
length for this line complex. The [Co III] A5893 arises from the 
3d 7 a 4 F- 3d 7 a 2 G multiplet, and is actually a blend of two lines — 
one at 5888.5 A and a second, but weaker, line at 5906.8 A (see Ap- 
pendix[A]and Table lAlt . Given the A values and wavelengths of the 
transitions contributing to the line complex, the weighted mean rest 
wavelength of the Co III line is 5892.7 A (note: this and previous 
are air wavelengths). Henceforth we use this value for calculating 
line velocities. 

3.2 Measuring the [Co ill] A5893 line flux 

For the main analyses in this work we focus on the flux in the 
[Co III] A5893 line. We measure the flux in this line as follows. 

We perform an initial Gaussian fit to the [Co III] A5893 line 
in order to the determine the center and width of the line. We then 
integrate the flux within ±1.5cr of the fitted line center and use 
this “integral” flux for the remainder of this paper. This integral 
boundary was chosen as a compromise between capturing a large 
fraction of the emitted line flux (97% for a strictly Gaussian profile) 
and limiting contamination from neighbouring emission lines. For 
SNe la with multiple nebular spectra we enforce common wave¬ 
length bounds for the flux integration at all epochs, as determined 
by the median fitted line center and width values across all epochs. 
Generally the integrated line flux and that calculated from the best 
Gaussian fit showed excellent agreement (see Figure[3}. but we pre¬ 
fer the integral flux as this is robust against non-Gaussianity of the 
line profile. 

To place our [Co III] A5893 line flux measurements on the 
correct absolute scale, we must ensure the spectra have the correct 
absolute flux calibration. To achieve this, we measure the expected 
B-band flux in the spectrum by convolving it with the B-band filter 
throughput curve and integrating. We then compute the ratio of this 
flux collected in the spectrum B passband to the true B-band flux 
of the SN at that epoch. The latter is determined from the late-time 
photometry for each of our SNe, as outlined in Appendix [Bl and 
presented in Table lB3l To ensure reproducability of our results, we 
report in this table the flux values derived from the raw measure¬ 
ments made from the spectra in their published form. 

We note that normalization with the B-band could introduce 
errors in the [Co III] A5893 flux due to chromatic errors in the spec¬ 
trum’s flux calibration. However, previous authors consistently per¬ 
formed chromatic flux calibration using spectrophotometric stan¬ 
dard stars, typically yielding excellent colour agreement with ob¬ 
served photometry (e.g. B~V scatter of 0.08 mag and 0.10 mag for 
the CfA and BSNIP samples, repectively). We also note that other 
systematic effects could affect our measurements of [Co III] A5893 
line flux. These include contamination from neighboring nebular 
emission lines (e.g. Fe II lines, see Figure O, residual host galaxy 
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Figure 2. Top panels: Comparison of radiative transfer (CMFGEN) model spectrum for the DDC10 llllondin et alj201~3l) delayed detonation model at very late 
epochs (left: +126 days, right: +300 days) compared to contemporaneous data for SN 201 lfe and SN 2012fr. Bottom panels: Emission spectra for various ions 
from CMFGEN for late-phase DDC10 models (epochs as above). 



Figure 3. Comparison of flux in the [Co III] A5893 line measured in two 
ways: strict integration of the spectrum flux within ±1.5cr of the fitted line 
center (y-axis and upper left inset), and the formal integral of the best fit 
Gaussian profile (x-axis and lower right inset). The solid line represents 
unity, while the dashed line is the mean ratio of the integral flux to Gaussian 
flux for the full sample (0.87 =t 0.05). 


light, or perhaps even p reviously undetected light echoes (see, e.g., 
ISpvromilio et alj20041 ). Thus we expect a conservative estimate for 
the systematic uncertainty in the [Co III] A5893 flux measurement 
to be about 10% of the measured flux. 

The final integrated [Co III] A5893 line flux, wavelength 
bounds for the integral, and synthetic B-band flux integrated from 
the spectrum are all presented in Table lB4l Variance spectra were 
not available for many of the literature SN la spectra in our analy¬ 
sis. To correct this, we smooth the spectrum with a Savitszky-Golay 
filter dSavitzkv & Golavll 19641) . then smooth the squared residuals 
of the data from this smooth curve to derive a variance spectrum 
mea sured direc tly from the noise in the data (as we did for data in 
IChildress et al.1120143 ). [Co III] A5893 line flux errors were then 
determined from these corrected variance spectra. 


4 EVOLUTION OF THE [Co ill] A5893 LINE FLUX 


4.1 Theoretical expectations for [Co ill] A5893 evolution 

The decay of 56 Co to o6 Fe produces positrons and energetic 
gamma-rays. The charged positrons carry kinetic energy which 
they lose to the surrounding medium via Coulomb interactions. 
At the nebular densities present at late times, the length scale 
for positron energy deposition is much smaller than the size of 
the nebula so the p ositrons essentially deposit all of their ki¬ 
netic energy locally dChan & Lingenfelteilll993l) . gamma-rays - 
either those emitted directly from 56 Co decay or created when 
the positrons annihilate - are subject to radiative transfer effects 
and will eventually free stream as the SN nebula expands and de¬ 
creases its density enough to become optically thin to gamma- 
rays. The onset of this phase - where positrons deposit a con¬ 
stant fraction of energy into the SN nebula and gamma-rays es¬ 
cape fully - has been observed i n late SN la bolometric light 
curves (e.g. Sollermanet alJl2004l: IStritzinger & Sollermanl l 2007l : 


iLeloudas et al.l2009l: Kerzendorf et alj2014bir 

Our expectation from a simple energetics perspective is that 
the flux of the [Co III] A5893 line should evolve as the square of the 
mass of cobalt as a function of time Mc 0 {t). The energy being de¬ 
posited into the nebula at these late phases arises from the positrons 
produced in o6 Co decay, and thus should scale with the mass of 
cobalt. If this energy is evenly deposited amongst all species in 
the nebula then the fraction of that energy absorbed by the cobalt 
atoms should be proportional to the mass fraction of cobalt. Thus 
the amount of energy absorbed by cobalt atoms follows the square 
of the cobalt mass as a function of time. If the fraction of that en¬ 
ergy emitted in the [Co III] A5893 line remains constant (see Sec¬ 
tion [6T2J then we expect a net quadratic dependence of the [Co III] 
A5893 line luminosity on the mass of cobalt as a funciton of time. 

Observational evidence for this temporal evolution of the 
[Co III] A5893 line should be expected from prior results. The late- 
phase bolometric light curves of SNe la closely follow the amount 
of ene rgy deposited by the decay of 5 6 Co (see, e.g., Scrilerman et al.l 


2004). It was also demonstrated by iKuchner et al 


19941) that the 


ratio of [Co III] A5893 to Fe 4700 emission follows the Co/Fe 
mass ratio (as noted above), and the Fe 4700 line flux generally 
scales with the total luminosity of the SN since Fe is the primary 
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coolant. These facts combine to lend an expectation that the net 
emission from the [Co III] A5893 line should scale quadratically 
with the mass of Co in the SN nebula as a function of time. Indeed, 
iMcClelland et ali d2013l) found such a quadratic dependence for the 
[Co III] A5893 line in SN 201 lfe. 

The above reasoning for Mq 0 dependence of the [Co III] 
A5893 flux holds for epochs when the nebula is fully transparent to 
gamma-rays. Thus it is important to inspect the theoretical expecta¬ 
tion for the timing of this gamma-ray transparency in the IGE zone. 
The energy released per decay from 56 Co is 3.525 MeV, of which 
3.3% is associated with the kinetic energy of the positrons, and we 
have ignored the energy associated with neutrinos. As the expan¬ 
sion is homologous, the optical depth associated with gamma-rays 
scales as 1/t 2 . Assuming that the kinetic energy of the positrons is 
captured locally, the energy absorbed per 56 Co decay in MeV is 

eco = 0.116 + 3.409 (l - exp[-T 0 (t 0 /t) 2 ]) (1) 

where r 0 is the effective optical depth at a time t 0 . If we denote f c as 
the time at which energy deposition by gamma-rays and positrons 
are equal, then EquationQ]can be rewritten as: 

Eco oc Mco (1 - 0.967 exp [-0.0346 (t c /t) 2 ] ) (2) 

We expect the flux from the [Co III] A5893 line would further scale 
as: 

p p _ xMco/Mige _ ... 

C C 1 + [a — 1 )Mco/Mige + bMot/MiGE 

where Mige is the total mass of the IGE zone, a and b are respec¬ 
tively the (time-dependent) factors relating the cooling efficiency 
of Co and other species (which have total mass of Mot) relative to 
iron, and x is the factor scaling the emission in the [Co III] A5893 
feature. If the thermal conditions in the SN nebula are relatively sta¬ 
ble (i.e. constant x) and cooling by non-iron species is negligible 
(i.e. the above denominator goes to unity), then the line flux simply 
becomes proportional to Mco/Mige ■ Combining Equations[2]and 
[3]yields: 

Fco oc (1 - 0.967 exp [-0.0346(f c /t) 2 ]) (4) 

For the DDC10 model, we find t c ~ 214 days (from explosion) 
- this would imply a deviation from Mo 0 of a factor of 2 from 
+ 150 to +400 days past maximum light (assuming a rise time of 
~ 17 days), or a factor of 1.5 from +200 to +400 days. Alternatively 
if t c ~ 80 days (see Section l4~3l > then the deviation from Mq 0 is 
only 20% from +150 to +400 days and 10% from +200 to +400 
days. 


4.2 Observed [Co ill] A5893 evolution in nebular spectral 
time series 

To examine the observed evolution of the [Co III] A5893 line, we 
turn to those SNe la with numerous nebular spectra. Specifically, 
we isolate the subset of SNe la in our sample with at least three 
epochs of observation later than +150 days past maximum. For the 
eight SNe la in our sample which meet this criterion, we also collect 
spectra between +100 < t < +150 days past maximum (dates 
listed in italics in Table Q}. These additional spectra allow us to 
further inspect the [Co III] A5893 flux evolution, but these spectra 
are not employed in our nickel mass estimates derived in Section[5] 
In the upper panel of Figure [4] we show the evolution of 
the [Co III] A5893 line luminosity versus time for our sample of 
SNe la with three or more observations after +150 days. We plot 
the line evolution for a linear (dotted) line and quadratic (solid 




Figure 4. Top: Evolution of the [Co III] A5893 line flux in SNe la with 
nebular time series (> 3 observations past 150 days), compared to curves 
following the mass of 56 Co as a function of time to the first power (dot¬ 
ted line) and second power (solid line). Data for each SN was shifted by 
a multiplicative offset (i.e. log additive offset) that reduced residuals with 
the Mco{t ) 2 line. Bottom: Evolution of the “Fe complex” flux with phase, 
compared to the same lines as above as well as an additional line propor¬ 
tional to the product of the 56 Co mass with the 5 ®Fe mass as a function of 
time (dashed curve). 


line) dependence on Mc 0 {t), with both curves normalized at phase 
t = +200 days. For each SN in this subset, we fit for a single multi¬ 
plicative scaling factor that minimizes the residuals of the Mcoif) 2 
line (i.e. we normalize each SN data set to that line - thus the rea¬ 
son for requiring multiple data points per SN). This isolates the 
time dependence of the line flux (which depends on the SN nebula 
physics) by removing its absolute magnitude (which depends on 
the quantity of 56 Ni produced). 

The evolution of the [Co III] A5893 line shows a remarkable 
agreement with the expected trend of Mc 0 {t) 2 , perhaps as early 
as phase +150 days. The one possible exception to the Mc 0 {t) 2 
trend is SN 199IT, which appears to have a shallower evolution 
than the other SNe la. As we show below (Section[+3j. this cannot 
arise from gamma-ray opacity. Instead the most likely explanation 
is probably a higher ionization state at early epochs (t < 300 days). 
Because of this, for SN 1991T only we excise epochs prior to 
300 days when calculating its 56 Ni mass in Section l+2l - a choice 
which yields more favorable agreement with previous analyses 
from the literature. 

To contrast the behavior of the [Co III] A5893 line with other 
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regions of the nebular spectra, we also inspected the evolution of 
the blue “Fe complex” of lines. For each spectrum we integrate the 
flux in the region 4100 — 5600 A (adjusted for each SN according to 
its central nebular velocity measured from the Co line) where the 
emission is almost entirely dominated by Fe lines (see Figure [2]>. 
Following our arguments for the expectation of the [Co III] A5893 
line flux, the Fe complex flux should be proportional to the energy 
being deposited - which scales as Mc 0 (t) - and the mass frac¬ 
tion of Fe (which should be relatively constant as Mco -C MF e 
at this point). Thus the Fe complex flux should scale linearly with 
Mco (t ). In the lower panel of Figure[4]we plot the evolution of the 
Fe flux for the sample of SNe la, and see that it follows more closely 
the Mco{t) curve than the Mc 0 (t ) 2 curve. However, we do note 
deviation from this line such that the logarithmic slope is some¬ 
what intermediate between 1 and 2. Additionally, earlier epochs are 
subject to a complicated interplay of additional energy deposition 
from gamma-rays (as for the [Co III] A5893 line, see Section l431 >. 
decreased emission due to nonzero optical depth in this region of 
the spectrum, and possible emission from Co II (see Figure^- 

We note that the above re sults a lso explain one aspect of the 
data presented in lForster et alJd2013l) . Those authors examined the 
late (35 < t < 80 days) colour evolution (i.e. Lira law) for a large 
sample of nearby SNe la and its relationship with dust absorption 
(as inferred from narrow sodium absorption). The mean value of B- 
band decline rates were roughly 0.015 mag/day, while the V-band 
decline rates were nearly twice that (0.030 mag/day). The B-band 
is dominated by the Fe complex whose flux decays as Mc 0 (t), 
while the V-band is heavily influenced by Co lines (see Figure[2]in 
Section im whose flux decays as Mco (t ) 2 . This naturally explains 
why the luminosity decay rate (in mag/day) in V-band is nearly 
twice that of the B-band, and contributes to why SNe la become 
bluer (in B — V) with time at these epochs. 


4.3 Testing gamma-ray opacity effects on [Co ill] A5893 
evolution 

While the data appear to agree with an Mc 0 {t) 2 dependence of 
the [Co III] A5893 flux evolution, it is important to investigate 
the impact of gamma-ray energy deposition on deviation from this 
parametrization. 

To this end, we isolated the subset of SNe la from our sample 
with at least one nebular spectrum earlier than +150 days and at 
least one spectrum later than +250 days. For the six SNe la satis¬ 
fying these criteria, we fit the [Co III] A5893 flux evolution using 
the parametrization of Equation|4] This fit has two free parameters: 
a multiplicative scaling for all the line fluxes, and the gamma-ray 
“crossing” time t c when energy deposition from gamma-rays and 
positrons are equal. These fits are shown in Figure[5] 

In general the [Co III] A5893 evolution is extremely well fit 
by this model, especially for SNe la with good temporal coverage 
and high signal-to-noise data (notably SN 201 lfe and SN 2012fr). 
Some SNe la have a gamma-ray crossing time similar to the pre¬ 
diction from our model ( t c ~ 200 days) while some other SNe la 
have shorter crossing times ( t c ~ 80 days). The implications of 
this for SN la progenitors will be discussed in further detail in Sec¬ 
tion [6T| Given these gamma-ray opacity model fit results, we cal¬ 
culate that deviations of [Co III] A5893 flux evolution from the 
simple Mco(t) 2 could range front 15% to 100% at / = 150, and 
7% to 55% at / = 200, and 4% to 30% at t = 250 days. 





100 150 200 250 300 350 400 

Phase (t) [days] 




Figure 5. Fits of gamma-ray opacity model to select SN la [Co III] A5893 
line fluxes. The fitted “crossing” time (when gamma-ray and positron en¬ 
ergy deposition are equal) is shown in each panel. 


5 MEASURING 56 NI MASS FROM SN la NEBULAR 
SPECTRA 

5.1 Placing [Co ill] A5893 flux measurements at disparate 
epochs on a common scale 

To place all our SN la [Co III] A5893 fluxes on a common scale, 
we first convert the observed line flux to the absolute line lumi¬ 
nosity emitted by the SN using the distance to the SN host galaxy. 
For some SNe la in our sample, redshift-independent distance mea¬ 
surements exist for the host galaxy, particularly a number with 
Cepheid distance measurements. For most of the SNe la in our 
sample, however, the SN distance is computed by converting the 
host galaxy redshift to a distance using a Hubble constant value 
of Ho = 73.8 km s _1 Mpc -1 chosen from lRiess et alj 1201 ll) to 
maintain consistency with those hosts with Cepheid distances from 
that work. For hosts with redshift-based distances, we assign a dis¬ 
tance uncertainty corresponding to a peculiar velocity uncertainty 
of 300 km s -1 . Table IbTI lists the full set of distance moduli (and 
references) employed in our sample. 

Calculating the absolute [Co III] A5893 flux emitted by each 
SN also requires correction for extinction by interstellar dust 
in the SN host galaxy. We accomplish this by calculating the 
ICardelli et al] d 19891 hereafter CCM) reddening curve at the rest 
central wavelength of the [Co III] A5893 complex for an appro¬ 
priate value of the reddening E(B — V ) and selective extinction 
Rv- For most SNe la in our sample, the reddening is extremely 
low (B(B — V) < 0 .10 m ag) , so w e use the light curve colour 
fitted by SiFTO ( IConlev et al. I2008I) . and a selective extinction 
value of Ry = 2.8 (appropriate for cosmological SNe la, see 
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Figure 6. Evolution of the absolute [Co III] A5893 line luminosity as a 
function of phase for all SNe la in our sample. The solid line corresponds 
to the square of the mass of 56 Co as a function of time, anchored by the 
[Co III] A5893 luminosity for SN 201 lfe at ~ 200 days. Here thick error 
bars correspond to flux measurement errors, while narrow error bars corre¬ 
spond to distance uncertainties. 


Figure 7. [Co III] A5893 line luminosity scaled to its equivalent value at 
t = 200 days using the Mp 0 (f) 2 curve (Tco) versus SN light curve 
stretch. As in Figure [6] thick error bars correspond to flux measurement 
errors, while narrow error bars corresponding to distance uncertainty aris¬ 
ing from peculiar velocities. 


IChotard et alfcOlU) . We note that the choice of Rv has negligible 
impact on the majority of our sample. SN la light curve colours 
are affected by both intrinsic colour and host galaxy extinction 
(see, e.g JScolnic et al.ll20l4) . so for SNe la with negative SiFTO 
colours - indicating blue intrinsic colours - we apply no colour 
correction (i.e. colour corrections never redden the data). In this 
work, we are not trying to standardize SN la (in which applying a 
colour correction to the intrinsic colours may also be appropriate); 
rather we are only concerned with eliminating the effects of dust 
extinction. 

Two SNe la in our sample, however, have strong extinction 
by unusual dust and thus must be treated differently. SN 2014J oc¬ 
curred behind a thick dust lane in the nearby starburst galaxy M82. 
[Foley et ail d2014l) performed a detailed fit of multi-colour photom¬ 
etry of the SN, and find it is best fit by a CCM-like reddening curve 
with E(B — V) = 1.19 and Rv = 1.64. We adopt their colour cor¬ 
rection for SN 2014J, and for the line flux uncertainty arising from 
the reddening correction we adopt their uncertainty for the visual 
extinction of ga v = 0.18 mag. SN 20071e showed moderately low 
extinction but with some variability in the sodium absorption fea¬ 
ture likely arising from interaction of the SN with its circumstellar 
medium I lSimon et al.ll2009l) . Despite this variability, most of the 
absorption strength remains stable, so we adopt a colour correction 
for SN 20071e with E (B - V) = 0.277 and R v = 2.56 as derived 
bv lsimon et al. MU). 

Figure[6]presents the total emitted [Co III] A5893 luminosity 
as a function of phase for all nebular spectra in our final sample. In 
this and subsequent figures, the thick errorbars represent the com¬ 
posite measurement errors from the [Co III] A5893 flux, B-band 
flux in the spectrum, observed (photometric) B-band magnitude, 
and extinction correction; the narrow error bars represent the dis¬ 
tance uncertainties. Points are colour-coded (in groups) based on 
the light curve stretch. 

The line luminosity values are then used to compute an effec¬ 
tive luminosity of the [Co III] A5893 line at a common phase of 
+200 days for all SNe la in the sample (hencefoward we refer to 
this as Leo) using the M^ 0 curve. For a single nebular spectrum, 


this can be calculated directly as: 

log(L Co (200)) = log(L Co (f)) + 7.80 x 10“ 3 * {t - 200) (5) 

For SNe la with multiple spectra, Leo is calculated as the \ 2 - 
weighted mean value across all acceptable epochs (150 < t < 
400 days) using the above equation. We note the above equa¬ 
tion is calculated assuming a time between explosion and B- 
band peak (i.e. rise time) of 17 days, but there may be an as- 
sociated uncertainty on this due to diversity in SN la rise times 
dGaneshalingam et al. 20 id) and p ossible dark phase before first 
light escapes dPiro & Nakaii i2013). Each day of difference in ex¬ 
plosion date results in a corresponding change in the final [Co III] 
A5893 luminosity of 1.8% - assuming an explosion date uncer¬ 
tainty of about 3 days, we thus expect the explosion date uncer¬ 
tainty contributes about 5% uncertainty to the final nickel mass de¬ 
rived in Section [5~2l 

As noted in Section[4] SN 199IT may represent a case where 
the stable ionization state is not established until later than other 
SNe (also evident in Figure [6j, so for this SN we use the later 
two epochs (t > 300) to establish Leo■ This also yields a favor¬ 
able agreement of our 56 Ni mass with literature estimates (see Sec¬ 
tion [5T2J. 

In Figure [7] we show the scaled t = 200 d [Co III] A5893 
line luminosity plotted against light curve stretch. A clear corre¬ 
lation is evident between the [Co III] A5893 line luminosity and 
stretch - this is expected given the [Co III] A5893 luminosity traces 
the amount of 56 Ni produced in the explosion, and o6 Ni directly 
powers the peak luminosity which correlates with the light curve 
stretch. 


5.2 Inferring Mjvi from [Co ill] A5893 flux 

Scaling the [Co III] A5893 flux values to the same phase (f = 
200 days) effectively places all measurements at the same epoch 
since explosion, so the amount of 56 Co will have the same propor¬ 
tionality to the amount of 56 Ni produced in the explosion. The final 
critical ingredient for inferring o6 Co mass (and thus 56 Ni mass) 
from the [Co III] A5893 line flux is the scaling between 56 Co mass 
and [Co III] A5893 flux. For reasons we will explore in Section [6)2l 
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we expect this conversion factor to be relatively stable in time (for 
phases 150 < t < 400 days considered here) and consistent across 
all SNe la. At these phases we also expect 56 Co to be the dominant 
isotope (by mass) of cobalt dSeitenzahl et al.ll2009bl) . a s 57 C o o nly 
dominates energy deposition around t ~ 1000 days dCiraur et~al1 
|20I3) . 

We expect the [Co III] A5893 line flux at phase t = 200 days 
to be linearly proportional to the mass of 56 Ni produced in ex¬ 
plosion (since the 56 Co mass fraction at this same epoch is nec¬ 
essarily the same for all SNe la). To convert [Co III] A5893 flux 
to S0 Ni mass requires some scaling between the two quantities to 
be determined. In principle this could be computed through radia¬ 
tive transfer modelling of late phases for SN la explosion mod¬ 
els. However, for simplicity in this work, we choose to anchor 
the relation with the well-studied SN la SN 201 lfe. Modelling of 
the ph otospheric phase light curve for SN 201 lfe bv IPereira et al.1 
d2013l) yielded a 5( ^Ni ma ss of M Ni = 0.53 ± 0.11 Mq. Re- 
centlv IMazzali et all (20151) extended their spectroscopic modelling 
of the SN 201 lfe spectral time series (presented for photospheric 
epochs in lMazzali et al.|[2014l ) to nebular phase epochs and found 
Mni = 0.47 ± 0.08ATq. For simplicity in this work, we thus will 
choose a sc Ni mass anchor for SN 201 lfe of Mjvi = 0.50Mq, 
yielding final 56 Ni mass values derived as: 

M Ni = 0.50 M 0 -^- (6) 

ill fe 


where log(Z/n/ e ) = 39.410 is the scaled [Co III] A5893 luminos¬ 
ity we measure for SN 201 lfe - this is used as a zeropoint for the 
remainder of our SN sample. The values for Mivt for our sample 
are presented in Table [3] In Section 1631 we further discuss the im¬ 
plications of our 5<J Ni mass values and their relation to the ejected 
masses of our SN la sample. 

Other techniques have been presente d for measuring the mass 
of 56 Ni produced in the SN la explosion. IStritzinger et al.l d2006ai) 
employed semi-empirical modelling of SN la bolometric light 
curves to measure the ejected mass and 56 Ni mass for a sample 
of 17 nearby SNe la. They then foun d that 56 Ni masses derived 
from modelling of the nebular spectra iMazzali et al. 1997, 


Stehle et aD l2005h yielded consistent results ( IStritzinger et al.1 


2006bl) . Seven of the SNe la from their sample are included in 


ours, and we show a comparison of our 5 °Ni values versus those 
derived from their two methods in Figure [8] In some of the cases, 
our 56 Ni masses are somewhat lower than theirs (both for the light 
curve and nebular 56 Ni mass estimates) though generally show ac¬ 
c eptable agreeme nt. We note that for SN 1994ae and SN 2002er 
. IStritzinger et ail <t2006bh employ a much higher reddening value 
than ours ( E(B — V) = 0.15 mag versus E(B — V) = 0.00 mag 
for SN 1994ae and E(B — V) = 0.36 mag versus E(B — V) = 
0.12 mag for SN 2002er), which is likely the source of the discrep¬ 
ancy between our values. 


6 DISCUSSION 

In this Section we discuss the important physical implications of 
our observational results above. First, we examine the fact that the 
[Co III] A5893 line flux evolution requires a constant scaling be¬ 
tween energy released by 56 Co decay and that absorbed by the 
nebula - this requires efficient local deposition of energy from 
positrons and near-complete escape of gamma-rays from the IGE 
core (Section f6.lt . Next, we argue that the [Co III] A5893 evolu¬ 
tion requires stable ionization conditions in the nebula for a period 


Table 3. Final SN la Nickel Masses 


SN 

ATjvi 



Mej 




(Mq) 

a 


(Mq) 

b 


SN1990N 

0.514 

± 

0.027(0.081) 

1.437 

± 

0.009 

SN1991T 

1.049 

± 

0.106(0.308) 

1.407 

± 

0.019 

SN1994ae 

0.458 

± 

0.013(0.069) 

1.417 

± 

0.013 

SN1995D 

0.593 

± 

0.059(0.165) 

1.448 

± 

0.009 

SN1998aq 

0.707 

± 

0.042(0.127) 

1.304 

± 

0.015 

SN1998bu 

0.686 

± 

0.029(0.292) 

1.299 

± 

0.027 

SN1999aa 

1.593 

± 

0.114(0.238) 

1.465 

± 

0.003 

SN2002cs 

0.775 

± 

0.081(0.130) 

1.361 

± 

0.016 

SN2002dj 

0.882 

± 

0.051(0.176) 

1.299 

± 

0.019 

SN2002er 

0.344 

± 

0.018(0.082) 

1.202 

± 

0.015 

SN2002fk 

0.625 

± 

0.016(0.120) 

1.346 

± 

0.016 

SN2003du 

0.414 

± 

0.022(0.177) 

1.373 

± 

0.010 

SN2003hv 

0.186 

± 

0.003(0.073) 

0.914 

± 

0.037 

SN2004bv 

1.294 

± 

0.040(0.266) 

1.468 

± 

0.003 

SN2004eo 

0.332 

± 

0.011(0.046) 

1.135 

± 

0.016 

SN2005cf 

0.625 

± 

0.044(0.184) 

1.308 

± 

0.013 

SN2007af 

0.440 

± 

0.029(0.071) 

1.289 

± 

0.017 

SN2007gi 

0.624 

± 

0.027(0.248) 

1.149 

± 

0.023 

SN20071e 

0.549 

± 

0.033(0.202) 

1.387 

± 

0.017 

SN2007sr 

0.609 

± 

0.027(0.107) 

1.311 

± 

0.045 

SN20091e 

0.673 

± 

0.065(0.102) 

1.380 

± 

0.026 

SN201 lby 

0.582 

± 

0.082(0.119) 

1.295 

± 

0.029 

SN201 lfe 

0.500 

± 

0.026(0.069) 

1.310 

± 

0.015 

SN201 liv 

0.349 

± 

0.046(0.122) 

0.818 

± 

0.032 

SN2012cg 

0.479 

± 

0.048(0.309) 

1.422 

± 

0.010 

SN2012fr 

0.670 

± 

0.043(0.287) 

1.454 

± 

0.004 

SN2012hr 

0.328 

± 

0.008(0.084) 

1.375 

± 

0.025 

SN2013aa 

1.658 

± 

0.091(0.717) 

1.468 

± 

0.004 

SN2013cs 

0.757 

± 

0.094(0.174) 

1.360 

± 

0.013 

SN2013dy 

0.608 

± 

0.047(0.137) 

1.450 

± 

0.004 

SN2013gy 

0.950 

± 

0.075(0.159) 

1.278 

± 

0.012 

SN2014J 

0.837 

± 

0.176(0.250) 

1.441 

± 

0.007 


“ Nominal uncertainties arise from measurement errors in the Co line flux 
or SN reddening, while distance uncertainties are listed in parenthesis. 
Systematic error for is estimated at 0.2Mg. 

b Includes only measurement uncertaintes from SN light curve stretch. 
Systematic error for M e j is estimated at 0.1 Mq. 


of several hundred days, which we support by demonstrating sta¬ 
bility of ionization-dependent flux ratios measured from the data 
(Section f6.2t . Finally, we discuss potential interpretations of SN la 
explosion conditions implied by our observed relationship between 
inferred j6 Ni mass and ejected mass (Section f6.3b . 


6.1 Gamma-ray transparency timescales for nebular SNe la 

For 56 Co to deposit a constant fraction of its decay energy into 
the nebula, positrons from the decay must be efficiently trapped 
in the IGE core and gamma-rays must be able to effectively es- 
ca pfl As noted above, efficient local positron energy deposition 
is expected to hold for t he temperatures and densities encountered 


at these nebular phases 

Axelrod 198(1 Chan & Lingenfelter 1993; 

Ruiz-Lapuente & Spruit 

1998). In practice, gamma-rays become 


3 We do note that other physical properties of the nebula (e.g. ionization 
or emission measure changes) could somehow conspire to compensate for 
gamma-ray opacity to make the line emission evolve as Mq 0 , but we con¬ 
sider the gamma-ray transparency scenario to be the simplest explanation. 
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Figure 8. Comparison of our 56 Ni mass values to those derived by 
iStritzinger et alj t 2006b ) using light curve modelling (left) and nebular 
spectra modelling (right). 


negligible after the time when the gamma-ray energy deposition 
equals that from positrons, which occurs when the optical depth 
drops enough to reach this equality (Section |4j}. We will refer to 
this henceforth as the “transparency” timescale t c - after this epoch 
positrons dominate energy deposition in the nebula. 

We fit the transparency timescale for several supemovae in 
Section 14.31 and found several have longer transparency times 
(t c ~ 180 days) close to the theoretical expectation for the 
DDC10 model. Previous analysis of gamma-ray transparency 
timescales found simil ar results: t c « 170 days for SN 2000cx 
dSollerman et alj 120041) and t c « 188 days for SN 2001el 
dStritzinger & Sollermanl 120071) . However we found that other 
SNe la (notably SN 201 lfe and SN 2007af) had much shorter 
transparency times ( t c ~ 80 days). This variation in transparency 
times may reflect a diversity in nebular densities, as most of the 
gamma-ray opacity at these late epochs will come from opacity 
from electrons in the nebula. Interestingly, the SNe la with shorter 
transparency times (SN 20life and SN 2007af) have lower stretch 
values than most of the SNe la with longer transparency times 
(SN 2003du, SN 2012fr, SN 2013aa), possibly indicating some re¬ 
lationship between nebular density and stretch. The one exception 
to this is SN 2003hv, which appears to have low stretch but long 
transparency time (and thus would im ply high density) - this re¬ 
sult is opposite to the findings of lMazzali et alj d201ll) who found 
SN 2003hv had reduced density in the inner regions of the ejecta. 
The source of this discrepancy is unclear, but may constitute further 
evidence that SN 2003hv is a “non-standard” event. 

Because of the diversity in gamma-ray transparency 
timescales in the SNe la we tested, it is likely that the impact of 
gamma-ray energy deposition on the [Co III] A5893 flux will be im¬ 
pacted by similar variability. Given the results above (Section [4.3t 
this may result in an average uncertainty of 30% on the final 56 Ni 
masses we infer. The only robust way to account for gamma-ray 
opacity effects is to obtain a nebular time series. However the trans¬ 
parency time is best constrained by observations from 100-150 days 
when the SN is only 3-4 magnitudes fainter than peak. Thus it 
should be observationally feasible to obtain such data for future 
SNe la observed in the nebular phase. 

More interestingly, the time evolution of the [Co III] A5893 
flux presents a new method for measuring the gamma-ray trans¬ 
parency time scale, as it gives a direct probe of the energy be¬ 
ing deposited into the nebula. Previously t his could only be done 
with the aid of bolometric light curves i Sollerman et alj 1 20041 : 


IStritzinger & Sollermanll2007l : lLeloudas et alJI 2009h . which neces 

sarily rely on extensive optical and infrared photometry and/or un¬ 
certain bolometric corrections. Instead, our method requires only 
two nebular spectra with contemporaneous optical photometry. 



Stretch 

Figure 10. Integrated flux ratios of the Fe III 4700 A complex compared 
to the Fe 11 5270 A complex (top) and Fe II 4200 A complex (bottom) as 
a function of light curve stretch (SiFTO s) for all SNe la in our sample. 
Formal spectrum flux error bars are smaller than the data markers. Markers 
are the same as for Figure[7] 


6.2 Ionization conditions in the SN nebula 

As noted above, the consistency of the [Co III] A5893 flux evolu¬ 
tion with the square of the cobalt mass implies a constant scaling 
between the energy being absorbed by cobalt atoms and the energy 
they emit in the [Co III] A5893 line. This implies stability in the 
ionization conditions of the nebula, which we now investigate from 
a more detailed inspection of our nebular spectra. 

To confirm that the ionization state of the nebula is indeed 
slowly evolving from phases 150 < t < 400 days, we examine the 
flux ratios of nebular emission lines arising primarily from Fe II and 
Fe III. If the ratio of these lines evolves with time, this would indi¬ 
cate a change in the ionization state. In the left panels of Figure [9] 
we highlight the regions of the typical SN la nebular spectra (here 
from SN 201 lfe and SN 2012fr) which are dominated by strong 
line complexes of either Fe II or Fe III. We integrate the flux in 
these regions for all the nebular SN la spectra in our sample, and 
in the right panels of Figure we show how the line flux ratios 
evolve with phase for the nebular time series SNe la (the same as 
from Section [4}- For this analysis we only consider phases later 
than t ~ 200 days, as this is when this region of the spectrum is 
reliably optically thin (see Section |4) - note this cuts SN 2007af 
from the Fe time series sample. 

Though there is indeed some evolution in the flux ratio of 
Fe II lines to Fe III lines, it is comparatively small - generally less 
than 10% change of the relative line flux in Fe III compared to 
Fe II. In sharp contrast, consider the Fe IIl/Fe II line flux ratios as 
measured from the t ~ 1 000 days spectrum for SN 201 lfe from 
iTaubenberger et alj J20I5I) - 0.52 for 4700/5270 versus a mean of 
1.6 at earlier phases, and 0.87 for 4700/4200 versus an earlier mean 
of 2.3 - which decrease by at least 65% from their values in the 
150 < t < 400 day range [we note these values should be consid¬ 
ered upper limits as it appears that the Fe III 4700 line has effec¬ 
tively disappeared in the t ~ 1000 days spectrum for SN 201 lfe, 
so the flux we measure here is likely due to other species]. By these 
very late phases the physical conditions in the SN la nebula have 
clearly changed in a dramatic fashion. Such is not the case for the 
SNe la in our sample at phases 150 < t < 400 days. 

In order to meaningfully compare the [Co III] A5893 line flux 
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Figure 9. Left panels: Multiple nebular phase spectra of SN 201 lfe (top) and SN 2012fr (bottom), highlighting the flux integration regions for the line 
complexes dominated by Fe II (light red regions) and Fe III (light magenta regions). Integration zones are shifted by the central redshift of the nebular Fe lines, 
here a blueshift of ~ 600 km s —1 for SN 201 lfe and a redshift of ~ 1800 km s —1 for SN 2012fr. Right panels: Temporal evolution of the ratios of the flux 
integral for the Fe 111 4700 A complex compared to the Fe II 5270 A complex (top) and Fe II4200 A complex (bottom) for SNe la with nebular time series. 


from different SNe la, another key requirement is that the ioniza¬ 
tion state of all SNe la be relatively similar. To test this assump¬ 
tion, we again use the Fe line ratios described above, but plot the 
mean Fe IIl/Fe II line flux ratio (computed as the error-weighted 
mean for SNe la with multiple epochs) versus light curve stretch in 
Figure [TO] We have excluded the highly reddened SN 20071e and 
SN 2014J to avoid any biases in these ratios due to uncertainty in 
the dust law (i.e. Rv). 

Flere we see some mild coherent change in the Fe line flux 
ratios (and thus ionization state) as a function of light curve stretch 
(with SN 2003hv as an outlier, as previously noted bv lMazzali et al .1 
1201 ll ). Flere the overall range of the line ratios is somewhat larger, 
with variations perhaps up to 40% but with a scatter of 7% (for 
4700/5270) and 15% (for 4700/4200). The ionization potentials of 
Fe and Co are very similar, which means a change in Fe III line flux 
induced by variation of the ionization state will manifest a compa¬ 
rable change in Co III line flux. Thus our [Co III] A5893 line fluxes 
above should have an additional scatter due to ionization state vari¬ 
ations of about 10%. Since our inferred 56 Ni masses are propor¬ 
tional to this line flux, this means that ionization state variations 
could induce a scatter of similar magnitude in our 56 Ni masses. 

Our measurement of the [Co III] A5893 line flux evolution, 
and variations of Fe IIl/Fe II line flux ratios as a function of both 
phase and SN stretch, coherently indicate that the ionization states 
of normal SNe la are remarkably consistent across different SNe 
and nearly constant across phases 150 < t < 400 days. This sta¬ 
bility of the ionization state was predicted bv lAxelrod d 1980t ), and 
our results here present the most compelling evidence to-date in 
support of that prediction. 


6.3 The relationship between o6 Ni and Ejected Mass 

The relationship between [Co III] A5893 luminosity and light curve 
stretch (Figure 0 hints at a relationship between physical proper¬ 
ties of the SN la progenitor system. In Section 15.21 we converted 
our measured [Co III] A5893 line luminosities into inferred o6 Ni 
masses. Here we convert light curve stretch into the SN ejected 
mass (i.e. progenitor mass for SNe la) using the relationship be¬ 
tween Me 2 _and_light curve stretch discovered by Scal zo et~aT1 
d2014al) . IScalzo et alJ d2014bt) used Bayesian inference to model 
the intrinsic distribution of ejected masses, which can be folded in 
as an additional prior when determining ejected mass using this re¬ 
lation. We derive a cubic fit to the relationship between stretch and 
M e j\ 

M ej = 2.07-7.51s + 11.56s 2 -4.77s 3 (7) 

= 1.35 + 1.30(s - 1) - 2.75(s - l) 2 - 4.77(s - l) 3 

The resultant values for ejected mass ( M e j ) we derive are presented 
in Table[3]along with our 5G Ni masses. 

In Figure [TT| we plot our inferred 56 Ni masses against these 
ejected masses. We note that there is a systematic uncertainty as¬ 
sociated with M e , calculation of about O.IMq, as determined by 
IScalzo et al.1 d2014al) from recovering masses of SN la explosion 
models. For 56 Ni masses, we previously noted several sources of 
uncertainty: 10% uncertainty in the [Co III] A5893 flux itself (Sec¬ 
tion 13.21 , 5% uncertainty on the t = 200 [Co III] A5893 lu¬ 
minosity due to uncertainty in the explosion date (Section [5Tb , 
10% from ionization state variations (Section E3, and possibly 
30% from variations in gamma-ray transparency timescales (Sec¬ 
tion [64}. Collectively this constitutes a possible 35% uncertainty 
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Figure 11. Mass of 56 Ni inferred from the scaled [Co III] A5893 line luminosity versus total ejected mass (i.e., progenitor WD mass) inferred from SN 
light curve stretch. As in Figure[6] thick error bars correspond to flux measurement errors (including colour correction uncertainties), while narrow error bars 
corresponding to distance uncertainties arising from peculiar velocities. The typical systematic uncertainties for estimating M e j (±0.1 M© - from the stretch- 
M e j relation of lScalzo~et all2014bh and Mffi (±O.2M 0 - from gamma-rays, explosion date uncertainty, and possible line contamination) are shown as the 
blue error bars in the left side of the plot. Note the anomalously high ATy, values for SN 1999aa and SN 2013aa, which we attribute to line contamination and 
distance uncertainty, respectively (see text for details). 


in our 56 Ni masses, which given the values we find would produce 
a mean uncertinty in Mjvt of about 0.2 Mq. 

The relation between Mni and M e j shows potential evidence 
for two regimes for the production of j6 Ni in SNe la. For sub- 
Chandrasekhar ejected masses (My < 1.3 M 0 - though note 
SN 1991bg-like objects are not included in this analysis), the 
amount of 56 Ni produced is clustered around Mm ~ O.4M 0 , with 
a possible increase of Mjvi with My (though we note the statistics 
are small). Chandrasekhar-mass progenitors (M e j « 1.4±O.1M 0 ) 
produce 56 Ni masses ranging from O.4M 0 < A/jvi < 1.2M 0 , 
with the extreme high 56 Ni masses (Mn% > 1-0M 0 ) occuring in 
SNe la spectroscopically similar to the peculiar SNe SN 1991T 
(SN 1999a a, SN 2004bv SN 2 013aa, and SN 199IT itself). 
Recently, iFisher & Jumped <1 20151) suggested that Chandrasekhar- 
mass SN la progenitors preferentially lack a vigorous deflagration 
phase following the initial ignition, and result in a nearly pure det¬ 
onation that produces about l.OM 0 of o6 Ni and shows similarity 
to SN 1991T. Our findings that the [Co III] A5893 luminosity is ex¬ 
ceptionally high only in 91T-like SNe la could lend support to this 
theory. 


We note that SN 1999aa and SN 2013aa have anomalously 
high Mjvi values (indeed exceeding their My values). We visually 
inspected the spectra of these SNe, and find no fault in our fits to the 
[Co III] A5893 line. SN 1999aa notably has the broadest linewidth 
of our sample, which could result in contamination of our measured 
[Co III] A5893 flux by nearby Fe II lines (see Figure[2]i. SN 2013aa 
has a relatively uncertain distance to its host galaxy. We expect the 
true Mnz for these two SNe is likely to be closer to that of the other 
SN 1991T-like SNe la, near 1.0-1.2M 0 . 

To compare model predictions with our inferred o6 Ni mass 
values, we gather ejected mass and 56 Ni mass outcomes from nu¬ 
merous SN la explosion models and plot them against our data in 
Figure |T2] These models can be generally grouped into three cate¬ 
gories: sub-Chandrasekhar mass detonations, Chandrasekhar-mass 
deflagration to detonation transitions (DDT), and Chandrasekhar- 
mass deflagrations which fail to detonate. We discuss each category 
and its agreement with the data below. 

Sub-Chandrasekhar (sub-Ch) mass detonations: We con¬ 
sider sub-Ch detonations from lsim et alj|201ol) . where detonations 
were artificially initiated in WDs of varying initial masses. These 
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Figure 12. Mjji and M e j values inferred from data (dark grey diamonds - 
data as in Figure II It compared to various the oretical mode ls: pure detona¬ 
tion models from sub-Chandra detonations (Sim et a lj20ld yellow circles) 
and double-detonations iFink et al]l201(l magenta d iamonds); detonation- 
defla gration transi t ions ( DDT) from I Seitenzahl et alj l2013bl blue sqaures) 
and lBlondin et alj d2013L cyan downward triangles); and Chadra-mass pure 
deflagrations IFink et alj2014l green upward triangls). The observed width- 
luminosity re l ation (and its scatter) from the recent cosmology analysis of 
[Betoule et aljt2014h are shown as the red curve (and light red shaded area). 


models are also applicable to sub-Ch WDs ignited via other mech¬ 
anisms (e.g. a violent merger), and were also employed to estimate 
the brightness distribution of violent mergers in iRuiter et al. 1 20131). 


We also examine sub-Ch double detonation mod els fr om 
(2010): these are qualitatively similar to the ISim et al.l 


Fink et alj 

20 id) but 


the ignition mechanism naturally arises from a surface helium layer 
ignition. Both models show a similar relationship between Mjvi 
and M e j, which shows a much steeper increase of Mn% with M e j 
than we infer from our data. However we note that with the sys¬ 
tematic uncertainty in M e j estimates (from stretch) these may be 
compatible with the data. 

Deflagration to detonation transitions (DDT): We present 
model s from both [Seitenzahl et al. (201 3bl) and iBlondin et al.l 
J20 1 3l) - including the DDC10 model employed for radiative trans¬ 
fer calculations in Section lXTl In general for these models the Mch 
progenitors undergo an initial deflagration phase which transitions 
to a detonation at a later time: the timing of this transition directly 
sets the amount of 56 Ni produced. For lSeitenzahl et al] d2013bl) . the 
DDT time was cal culated from the sub-grid scale turbulent energy 
ICiaraldi-Schoolmann et alj2013l ) which in practice varied with the 
vigorousness of the initial defla grat ion (set by han d as the number 
of initial ignition points). For lBlondin et al.l <12013h . the DDT time 
is set by a manual trigger. Both sets of models cover a range of 56 Ni 
mass production, similar to the range inferred from our data. 

Pure deflagrations: Finally we consider pure deflagration 
models presented in iFink et all 


20141). These models 


are varia¬ 


tions on the Chandrasekhar-mass I Seitenzahl et alj d2013bl) models 
in which the DDT module has been intentionally turned off. Many 
of these deflagration models fail to fully unbind the star and eject 
only a portion of the WD’s t otal mas s and leave a bound remnant - 
we note that the lScalzo et al .1 J2014bl) method for estimating ejected 
mass from light curves is not trained to account for bound remants 
so may have some additional systematic uncertainty for this explo¬ 
sion mechanism. Interestingly, these models show a weak depen¬ 


dence of Mjvi on Mej for sub-Ch ejected masses, similar to what 
we infer for this regime of the data. This also shows some agree¬ 
ment with the width-luminosity relation observed in cosmological 
supernova samp le s - w e show this as well in Figure [12] using the 
WLR from Beto ule et al.l 1 20j4 ) converted to M^j and M e j using 
the relations presented in lScalzo et al.l d2014bh . 

If indeed the M^i-Mej trend arises from two distinct ex¬ 
plosion mechanisms for SNe la, several key questions remain to 
be answered with future research. One such question is where the 
split between the two mechanisms occurs - SNe la at Mch with 
~ 0.5 Mq of 56 Ni could arise from either mechanism - and what 
physical property of the progenitor decides which mechanism oc¬ 
curs. Next we should investigate why the two mechanisms produce 
SNe la which obey the same width-luminosity relation - one might 
expect that a different relationship between Mjvi and M e j would 
yield different relationship between peak luminosity and light curve 
width. Such insights could be further advanced by study of other re¬ 
lated thermon uclea r explosions which span a broader range of Mni 
and M e j (e.g. jMcCullv et al1l2014bl . see their Figure 15). 

Finally, we most critically should assess whether the two 
mechanisms calibrate cosmological distances in the same fash¬ 
ion. Recent evidence has been mounting that SNe la show 
progenitor signatures (e.g. CSM interaction, high-velocity fea¬ 
tures, h ost galaxy properties ) which appear to clump into two 
groups iMaguire et alj 120141 ). In parallel, SN la cosmological 
analyses have found that SNe la in high- and low- mass g alax- 


2O10j; Kelly et alj 2010 

; Lanmeitl et al. 2010: Gupta et al. 

2011]; 

D'Andrea et al. 2011; 

Konishi et al. 2011 

; Galbanv et al. 

2012; 

Hayden et al. 2013: Johansson et al. 2013; 

Childress et al. 

2013 b; 

Rigault et al. 2013; Childress et al. 2014c; Kelly et al. 2015|) 

These 


and the current study motivate further examination of the environ¬ 
ments and standardized luminosities of SNe la whose 56 Ni mass 
and ejected mass are assessed with the techniques presented here. 
Such a study is limited by distance uncertainties, and thus should be 
targeted at SNe la in the nearby smooth Hubble flow (z > 0.015) 
where distance uncertainties from peculiar velocities become small 
(< 0.10 mag). 


7 CONCLUSIONS 

In this work we examine the [Co III] A5893 feature in 94 nebu¬ 
lar phase (150 < t < 400 days past peak brightness) spectra of 32 
SNe la compiled from the literature and new observations. This fea¬ 
ture arises predominantly from radioactive 56 Co, the decay product 
of 56 Ni (which powers the bright early light curve) - thus this fea¬ 
ture provides a direct window for investigating the power source 
behind SN la light curves. 

We used nebular time series for eight SNe la to show that the 
temporal evolution of the [Co III] A5893 flux falls very close to the 
square of the mass of 56 Co as a function of time. This is the ex¬ 
pected dependence in the limit where the nebula is fully optically 
thin to gamma-rays produced in the 56 Co decay but locally ther- 
malizes energy from positrons emitted in the decay. We then used 
this uniform time dependence to infer the relative amount of 56 Ni 
produced by all 32 SNe la in our sample by using SN 201 lfe as an 
anchor (at Mjvi = 0.5 Mq). 

The greatest systematic uncertainty in our 56 Ni mass measure¬ 
ments was the time at which the nebula becomes effectively opti¬ 
cally thin to gamma-rays (which we define by the “crossing” time 
when energy deposition from positrons begins to exceed that of 
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gamma-rays). Though this could intoduce 30% uncertainty in 56 Ni 
masses (on average, though this is time dependent), we showed that 
the gamma-ray transparency time can be readily measured when 
multiple nebular spectra are available. In particular, a single spec¬ 
trum at phases 100 < t < 150 days past maximum light - when the 
SN is only 3-4 magnitudes fainter than peak - can easily constrain 
the gamma-ray transparency time. This can robustify our technique 
for measuring 56 Ni masses of future SNe la, but the gamma-ray 
transparency time itself could provide important clues to SN la pro¬ 
genitor properties. 

When comparing our inferred 56 Ni masses to the ejected 
masse s of our SN la sample (using techniques from iScalzo et al.l 
l2014afi ). we find evidence for two regimes in the production of 
56 Ni (which are too distinct to be an artefact of systematic un¬ 
certainties). For low ejected masses (low stretch), Mni clusters 
at low values (Mjv; se 0.4 Mq). At high ejected masses (high 
stretch) near the Chandrasekhar mass, Mn! has a much larger 
spread (0.4 Mq < Mjvi < 1.2 Mq). This could constitute evi¬ 
dence for two distinct explosion mechanisms in SNe la. 

This work has illustrated the power of the nebular [Co III] 
A5893 feature in probing the fundamental explosion physics of 
SNe la. We provide a simple recipe for calculating S0 Ni mass from 
the [Co III] A5893 line flux from a single nebular epoch, as well as 
prescription for a more robust measurement that accounts for opac¬ 
ity effects by using multiple nebular epochs. Future measurements 
which eliminate the opacity systematics and distance uncertainties 
could provide a detailed understanding of the explosion mecha¬ 
nisms for SNe la. 
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Table Al. Atomic data for Co III 


Lower level 

Upper level 

A° 

A 

E(eV) 

3d 7 a 4 F 9/2 

3d 7 a 2 G7/2 

0.014 

5627.104 

2.203 

3d 7 a 4 F 9/ 2 

3d 7 a 2 G 9/2 

0.400 

5888.482 

2.105 

3d 7 a 4 F 7/2 

3d 7 a 2 G 7/2 

0.150 

5906.783 

2.203 

3d 7 a 4 F 7/2 

3d 7 a 2 G 9/2 

0.120 

6195.455 

2.105 

3d 7 a 4 F 5/2 

3d 7 a 2 G 7/2 

0.110 

6127.670 

2.203 


° The A values are from Hansen et_al.| j 19841) 



Figure Al. Illustration of the variation of the flux in the [Co III] A5893 
blend with electron density. The different curves are for temperatures of 
3000 K (red curve) to 10000 K (top curve) in steps of 1000K. For the cal¬ 
culation a SN age of 100 days and a N(CoIII) to N(Ne) ratio of 0.01 was 
assumed. The emission measure has been divided by the product of the 
electron and Co densities. 


APPENDIX A: ATOMIC DATA FOR Co ill 

Rest wavelength: Table IaTI lists the transitions in the multiplet of 
Co III contributing to the 5893 A feature. Note that the second and 
third transitions contribute to the main 5893 A, while other transi¬ 
tions in the multiplet produce other features of interest such as the 
~ 6150 A feature of Co III (see Figure[2}. 

Line emissivity: The [Co III] A5893 blend accounts for about 
70% of the multiplet cooling. Above an electron density of 10' 
cm' 3 and a temperature of 4000 K, the above multiplet accounts 
for between 50 and 70% of all cooling by Co III. As the tempera¬ 
ture and/or electron density is lowered below these values, cooling 
via the 3d 7 a 4 F- 3d 7 a 4 P multiplet (with an excitation energy of 
~1.9eV) becomes increasingly important. Plots of the flux in the 
[Co III] A5893 blend as a function of electron density and temper¬ 
ature are shown in Fig. IA1I The critical density for the 3d 7 a 2 G 
levels is ~ 3 x 10 7 and for the 3d' a 4 P levels it is ~ lx 10 6 . 
We note here that collision rates for Co II I are unavailable and so 
we adopted collision strengths from Sunde rland et al] J2002) which 
were computed for Ni IV, and which has a similar electronic struc¬ 
ture to Co III. 
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APPENDIX B: ANCILLARY DATA ON SN SAMPLE 


This Appendix presents the properties of SNe la and their host 
galaxies derived in the literature or measured from the data for our 
analysis. 

Table IB ll presents the host galaxy distance and redshift infor¬ 
mation for the SNe la in our sample, along with references. As 


Riess et al. 

(2011) 

Riess et al. 

(2011) 


erentially derived from Cepheid distances from 1 
where available, and from host redshift with the 1 
Hubble constant value (Ho = 73.8 km s -1 Mpc -1 ). Generally 
we avoid other distance indicators such as Tully-Fisher, surface 
brightness fluctuations, and planetary nebula luminosity functions. 
The_one exc eption to this is SN 2012cg, for which we use the 


Cortes et al. 000^1 Tully-Fisher distance modulus, as was done in 


Silverman et all J 2012bl) . 


For some SNe la (SN 2003du, SN 201 lby, SN 2013dy), pre¬ 
vious analyses in the literature have found that the redshift-based 
distance modulus yields anomalous values for the SN peak lu¬ 
minosity. Thus for these we e mploy SN-based distanc e moduli 
adopted by pre vious authors: t hefstritzinger et al. (2006b) value for 
SN 2003du, the Graham et alJ d2015h value for SN 201 lby, and the 
IPan et al.1 J20 1 5al) value for SN 2013dy (derived with SNooPy). 

Table lB2| presents the pertinent light curve fit parameters from 
the main ph otospheric p hase light curves. Light curves were fit 
with SiFTO dConlev et alJ 12008) . For some SNe la (SN 201 liv, 
SN 201 lby, SN 2012ht), photometry was not explicitly published, 
so we utilize the published light curve widths and colours for these 
SNe. Where appropriate, we convert Ami 5 (B) or SALT2 xl to 
SiFTO stretch s using the relations published in lGuv et alJ d2007l . 
120101 ) . 


Light curves for SN 2013cs and SN 2013gy were obtained 
with th_e Las Cumbres Observatories Global Telescope (LCOGT 
iBrown et all2013l ) and have been reduced using a custom pipeline 
developed by S. Valenti. The pipeline employs standard procedures 
(PyRAF, DAOPHOT) in a Python framework. Host galaxy flux was 
removed using a low order polynomial background. Point spread 
function magnit udes were transfo rmed to the standard Sloa n Digi¬ 
tal Sky Survey dSmith et al .11 20021) filter system (for gri) or lLandoltl 
( 1992 1 system (for BV) via standard star observations taken during 
clear nights. These light curves will be released in a future LCOGT 
supernova program publication. 

Table lB3l presents the late-time photometry for our SNe la. 
The B-band magnitudes for each phase of spectroscopic observa¬ 
tion was determined from a linear fit to the very late (t > 60 days 
past peak) B-band light curve. Spectra front T15b were already flux 
calibrated to contemporaneous photometry, so here we report the 
B-band flux synthesized from the spectrum. Similarly, M15 spec¬ 
tra were flux calibrated with spectrophotometric standard stars, and 
thus are flux calibrated modulo any changes in grey extinction or 
slit loss between the standard star and SN observation. 

For some SNe la, sufficient late time photometry was not 
available to reliably extrapolate the B-band luminosity at the time 
of spectroscopic observations - generally this occurred if the SN 
did not have photometry beyond 100 days past peak or if a public 
light curve was not available. For these SNe, late photometry from 
an alternate SN was employed, with the requirement that the light 
curve stretch be exceptionally close. These SNe with their surro¬ 
gate light curves are shown in Figure IB 11 and the surrogate SN is 
denoted in the reference column of Table lB3l 

Table lB4l presents the measured [Co III] A5893 line fluxes, 
synthesized B-band fluxes, and wavelength integration bounds for 


the [Co III] A5893 flux measurements. Note that the fluxes in this 
table are in the natural units of the spectroscopic data, as measured 
from its publicly available format. 
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Table Bl. Nebular SN la Host Information 


SN 

Host 

Name 

Host z 

Host 
(mag) a 

Method 

z Ref. 

li Ref. 

SN1990N 

NGC 4639 

0.003395 

31.67 ±0.07 

Cepheids 

Wong et al. (2006) 

Riess et al. (2011) 

SN1991T 

NGC 4527 

0.005791 

30.72 ±0.13 

Cepheids 

Strauss et al. (1992) 

NED 

SN1994D 

NGC 4526 

0.002058 

29.73 ±1.06 

Redshift 

Cappellari et al. (2011) 


SN1994ae 

NGC 3370 

0.004266 

32.13 ±0.07 

Cepheids 

Krumm & Salpeter (1980) 

Riess et al. (2011) 

SN1995D 

NGC 2962 

0.006561 

32.25 ±0.33 

Redshift 

Cappellari et al. (2011) 


SN1998aq 

NGC 3982 

0.003699 

31.70 ±0.08 

Cepheids 

de Vaucouleurs et al. (1991) 

Riess et al. (2011) 

SN1998bu 

NGC 3368 

0.002992 

30.20 ±0.20 

Cepheids 

de Vaucouleurs et al. (1991) 

NED 

SN1999aa 

NGC 2595 

0.014443 

33.98 ±0.15 

Redshift 

de Vaucouleurs et al. (1991) 


SN2002cs 

NGC 6702 

0.015771 

34.17 ±0.14 

Redshift 

Trager et al. (2000) 


SN2002dj 

NGC 5018 

0.009393 

33.04 ± 0.23 

Redshift 

Rothberg & Joseph (2006) 


SN2002er 

UGC 10743 

0.008569 

32.84 ±0.25 

Redshift 

de Vaucouleurs et al. (1991) 


SN2002fk 

NGC 1309 

0.007125 

32.59 ± 0.09 

Cepheids 

Koribalski et al. (2004) 

Riess et al. (2011) 

SN2003du 

UGC 9391 

0.006384 

32.75 ±0.20 

SN 

Schneider et al. (1992) 

Stritzinger et al. (2006b) 

SN2003hv 

NGC 1201 

0.005624 

31.92 ±0.39 

Redshift 

Ogando et al. (2008) 


SN2004bv 

NGC 6907 

0.010614 

33.31 ±0.20 

Redshift 

Scarano et al. (2008) 


SN2004eo 

NGC 6928 

0.015701 

34.16 ±0.14 

Redshift 

Theureau et al. (1998) 


SN2005am 

NGC 2811 

0.007899 

32.66 ±0.28 

Redshift 

Theureau et al. (1998) 


SN2005cf 

MCG-01-39-3 

0.006430 

32.21 ±0.34 

Redshift 

Childress et al. (2013a) 


SN2007af 

NGC 5584 

0.005464 

31.72 ±0.07 

Cepheids 

Koribalski et al. (2004) 

Riess et al. (2011) 

SN2007gi 

NGC 4036 

0.004620 

31.49 ±0.47 

Redshift 

Cappellari et al. (20H) 


SN20071e 

NGC 7721 

0.006721 

32.31 ±0.32 

Redshift 

Koribalski et al. (2004) 


SN2007sr 

NGC 4038 

0.005477 

31.66 ±0.08 

Cepheids 

Lauberts & Valentijn (1989) 

Riess et al. (2011) 

SN2008Q 

NGC 524 

0.008016 

32.69 ±0.27 

Redshift 

Cappellari et al. (2011) 


SN2009ig 

NGC 1015 

0.008770 

32.89 ±0.25 

Redshift 

Wong et al. (2006) 


SN20091e 

MCG-04-06-9 

0.017786 

34.44 ±0.12 

Redshift 

Theureau et al. (1998) 


SN201 lby 

NGC 3972 

0.002843 

32.01 ±0.07 

SN 

Verheiien & Sancisi (2001) 

Graham et al. (2015) 

SN201 lek 

NGC 918 

0.005027 

31.67 ±0.43 

Redshift 

Rhee & van Albada (1996) 


SN201 lfe 

NGC 5457 

0.000804 

29.05 ±0.06 

Cepheids 

de Vaucouleurs et al. (1991) 

Shappee & Stanek (2011) 

SN201 liv 

NGC 1404 

0.006494 

32.23 ±0.33 

Redshift 

Graham et al. (1998) 


SN2012cg 

NGC 4424 

0.001458 

30.90 ±0.30 

TF 

Kent et al. (2008) 

Cortes et al. (2008) 

SN2012fr 

NGC 1365 

0.005457 

31.31 ±0.20 

Cepheids 

Bureau et al. (1996) 

Silbermann et al. (1999) 

SN2012hr 

ESO 121-G026 

0.007562 

32.56 ±0.29 

Redshift 

Koribalski et al. (2004) 


SN2012ht 

NGC 3447 

0.003559 

30.92 ±0.61 

Redshift 

Guthrie & Napier (1996) 


SN2013aa 

NGC 5643 

0.003999 

31.18 ±0.54 

Redshift 

Koribalski et al. (2004) 


SN2013cs 

ESO 576-G017 

0.009243 

33.00 ±0.24 

Redshift 

Pisano et al. (2011) 


SN2013dy 

NGC 7250 

0.003889 

31.49 ±0.10 

SN 

Schneider et al. (1992) 

Pan et al. (2015a) 

SN2013gy 

NGC 1418 

0.014023 

33.92 ±0.15 

Redshift 

Catinella et al. (2005) 


SN2014J 

NGC 3034 

0.000677 

27.60 ±0.10 

Ceuheids+SN de Vaucouleurs et al. (1991) 

Folev et al. (2014) 


“ For hosts with redshift-based //., uncertainty includes a peculiar velocity term of 300 km s 1 . 
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Table B2. Light Curve Fit Results and Sources 


SN 

MJD 

of Bmax 

Rest Frame 

Bmax (mag) 

Stretch 

Color c 

LC 

Ref. 

SN1990N 

48081.69 

12.706 ± 0.015 

1.081 ±0.012 

0.049 ±0.013 

Lira et al. (1998) 

SN1991T 

48373.95 

11.468 ± 0.028 

1.047 ±0.019 

0.101 ±0.026 

Liraet al. (1998) 

SN1994D 

49431.52 

11.775 ± 0.009 

0.823 ± 0.004 

-0.099 ±0.005 

Patat etal_£l 996) 

SN1994ae 

49684.56 

12.968 ± 0.019 

1.057 ±0.014 

-0.060 ±0.016 

CfA 

SN1995D 

49767.50 

13.273 ± 0.025 

1.097 ±0.014 

-0.006 ±0.015 

CfA 

SN1998aq 

50930.24 

12.316 ± 0.009 

0.965 ±0.010 

-0.146 ±0.009 

Riess et al. (2005) 

SN1998bu 

50952.06 

12.118 ± 0.013 

0.962 ±0.018 

0.270 ±0.010 

Jha et al. (1999); Suntzeff et al. (1999) 

SN1999aa 

51231.84 

14.755 ± 0.016 

1.134 ±0.009 

-0.047 ± 0.009 

CfA 

SN2002cs 

52409.24 

15.138 ± 0.039 

1.007 ±0.013 

0.017 ±0.016 

LOSS 

SN2002dj 

52450.32 

13.974 ± 0.035 

0.962 ±0.013 

0.098 ±0.016 

LOSS 

SN2002er 

52524.03 

14.267 ±0.057 

0.901 ±0.009 

0.123 ±0.018 

LOSS 

SN2002fk 

52547.28 

13.152 ± 0.017 

0.995 ±0.012 

-0.142 ±0.012 

LOSS 

SN2003du 

52765.48 

13.476 ± 0.007 

1.016 ±0.008 

-0.110 ±0.008 

LOSS 

SN2003hv 

52891.12 

12.444 ± 0.020 

0.741 ±0.021 

-0.115 ±0.014 

Leloudas et al. (2009) 

SN2004bv 

53159.83 

13.938 ± 0.024 

1.146 ±0.014 

0.122 ±0.013 

LOSS 

SN2004eo 

53277.66 

15.099 ± 0.038 

0.863 ± 0.009 

0.002 ±0.014 

CSP 

SN2005am 

53436.97 

13.698 ± 0.023 

0.710 ±0.019 

0.057 ±0.012 

CSP 

SN2005cf 

53533.28 

13.625 ± 0.035 

0.968 ± 0.009 

0.021 ±0.013 

LOSS 

SN2007af 

54173.88 

13.180 ± 0.015 

0.955 ±0.011 

0.058 ±0.010 

CSP 

SN2007gi 

54327.66 

13.158 ± 0.016 

0.871 ±0.013 

0.097 ±0.017 

Zhang et al. (2010) 

SN20071e 

54398.62 

13.876 ± 0.016 

1.028 ±0.015 

0.342 ±0.014 

CSP 

SN2007sr 

54448.70 

12.809 ± 0.042 

0.970 ±0.031 

0.141 ±0.016 

LOSS 

SN2008Q 

54505.47 

13.510 ± 0.033 

0.803 ± 0.025 

0.001 ±0.018 

LOSS 

SN2009ig 

55079.18 

13.437 ±0.019 

1.136 ±0.028 

0.055 ±0.019 

CfA 

SN20091e 

55165.16 

15.351 ± 0.030 

1.022 ±0.022 

0.079 ± 0.033 

CfA 

SN201 lby 

55690.60 

12.890 ± 0.030 

0.959 ± 0.019 

0.000 ± 0.050 

Graham et al. (2015) 

SN201 lek 

55788.90 

13.840 ± 0.150 

0.900 ± 0.020 

0.180 ±0.040 

Maguire et al. (2013) 

SN201 lfe 

55814.51 

9.940 ±0.010 

0.969 ±0.010 

-0.066 ±0.021 

Pereira et al. (2013) 

SN201 liv 

55906.00 

12.530 ± 0.040 

0.684 ± 0.020 

0.000 ± 0.050 

Foley et al. (2012) 

SN2012cg 

56082.03 

12.128 ± 0.011 

1.063 ±0.011 

0.184 ±0.010 

Munari et al. (2013) 

SN2012fr 

56243.68 

12.017 ±0.013 

1.108 ±0.008 

0.059 ±0.011 

Zhang et al. (2014) 

SN2012hr 

56289.20 

13.780 ± 0.020 

1.018 ±0.021 

0.030 ±0.010 

Maguire et al. (2013) 

SN2012ht 

56295.60 

13.220 ± 0.040 

0.818 ±0.025 

0.000 ± 0.050 

Yamanaka et al. (2014) 

SN2013aa 

56344.00 

11.330 ± 0.050 

1.146 ±0.019 

-0.050 ±0.010 

Maguire et al. (2013) 

SN2013cs 

56436.95 

13.659 ± 0.032 

1.006 ±0.010 

0.054 ±0.015 

LCOGT 

SN2013dy 

56500.65 

12.824 ± 0.048 

1.101 ±0.007 

0.195 ±0.015 

Pan et al. (2015b) 

SN2013gy 

56648.61 

14.742 ± 0.020 

0.948 ± 0.008 

0.070 ±0.011 

LCOGT 

SN2014J 

56689.74 

11.850 ± 0.012 

1.086 ±0.010 

1.251 ±0.012 

Folev et al. (2014) 


CfA: lRiess et alj jl999t) : fjha et alj j2006t) : [Hicken et al] j2009l[2Q12l) ; LO SS:lGaneshalingam et al] J2010I) ; CSP: IContreras et~afl j2O10l) ; IStritzinger et alj 

I201 Hi 
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Figure Bl. Surrogate light curves for SNe la with poor late time light curve coverage. Original SN data are shown as small black dots while surrogate SN data 
are shown as larger blue circles. 
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Table B3. Late Phase Photometry 


SN 


Phase t 
(days) 


Obs. Date 


mg (t) 
(mag) 


Phot. Ref. 


17.394 ± 0.014 Lira et at. (19981 


SN1990N 

160 

19901217 

SN1990N 

186 

19910112 

SN1990N 

227 

19910222 

SN1990N 

255 

19910322 

SN1990N 

280 

19910416 

SN1990N 

333 

19910608 

SN1991T 

113 

19910819 

SN1991T 

186 

19911031 

SN1991T 

258 

19920111 

SN1991T 

320 

19920313 

SN1991T 

349 

19920411 

SN1993Z 

173 

19940317 

SN1993Z 

205 

19940418 

SN1994D 

613 

19951124 

SN1994ae 

144 

19950422 

SN1994ae 

153 

19950501 

SN1995D 

277 

19951124 

SN1995D 

285 

19951202 

SN1996X 

298 

19970210 

SN1998aq 

211 

19981124 

SN1998aq 

231 

19981214 

SN1998aq 

241 

19981224 

SN1998bu 

179 

19981114 

SN1998bu 

190 

19981125 

SN1998bu 

208 

19981213 

SN1998bu 

217 

19981222 

SN1998bu 

236 

19990110 

SN1998bu 

243 

19990117 

SN1998bu 

280 

19990223 

SN1998bu 

329 

19990413 

SN1998bu 

340 

19990424 

SN1999aa 

256 

19991109 

SN1999aa 

282 

19991205 

SN2002cs 

174 

20021106 

SN2002dj 

222 

20030201 

SN2002dj 

275 

20030326 

SN2002er 

216 

20030410 

SN2002fk 

150 

20030227 

SN2003du 

109 

20030823 

SN2003du 

138 

20030921 

SN2003du 

139 

20030922 

SN2003du 

142 

20030925 

SN2003du 

209 

20031201 

SN2003du 

221 

20031213 

SN2003du 

272 

20040202 

SN2003du 

377 

20040517 

SN2003gs 

170 

20040214 

SN2003hv 

113 

20031228 

SN2003hv 

145 

20040129 

SN2003hv 

323 

20040725 

SN2004bv 

171 

20041114 

SN2004eo 

228 

20050516 

SN2005cf 

319 

20060427 

SN2006X 

127 

20060626 

SN2006X 

152 

20060721 

SN2006X 

277 

20061123 

SN2006X 

307 

20061223 


17.760 

± 

0.011 


18.337 

± 

0.007 


18.731 

± 

0.008 


19.083 

± 

0.011 


19.829 

± 

0.018 


17.026 

± 

0.032 

Schmidt et al. (1994) 

17.614 

± 

0.024 


18.194 

± 

0.018 


18.694 

± 

0.018 


18.928 

± 

0.019 


19.064 

± 

0.063 

Ho et al. (2001) 

19.588 

± 

0.094 


23.227 

± 

0.040 

Altavilla et al. (2004) 

17.672 

± 

0.025 

Altavilla et al. (2004) 

17.812 

± 

0.027 


19.887 

± 

0.041 

Altavilla et al. (2004) 

20.001 

± 

0.043 


20.108 

± 

0.022 

Salvo etal. (2001) 

18.140 

± 

0.011 

Riess et al. (2005) 

18.442 

± 

0.013 


18.593 

± 

0.013 


17.202 

± 

0.017 

SN2011fe 

17.347 

± 

0.019 


17.586 

± 

0.022 


17.705 

± 

0.023 


17.956 

± 

0.026 


18.049 

± 

0.028 


18.538 

± 

0.034 


19.187 

± 

0.042 


19.332 

± 

0.043 


20.872 

± 

0.000 

SN1991T 

21.236 

± 

0.000 


20.264 

± 

0.080 

LOSS 

20.143 

± 

0.058 

Pignata et al. (2004) 

21.003 

± 

0.040 


20.429 

± 

0.022 

Pignata et al. (2008) 

18.099 

± 

0.024 

CfA, LOSS 

17.612 

± 

0.004 

Stanishev et al. (2007) 

18.043 

± 

0.003 


18.058 

± 

0.003 


18.102 

± 

0.003 


19.096 

± 

0.004 


19.274 

± 

0.004 


20.031 

± 

0.005 


21.589 

± 

0.008 


19.877 

± 

0.051 

LOSS 

16.927 

± 

0.003 

Leloudas et al. (2009) 

17.436 

± 

0.004 


20.870 

± 

0.012 


18.713 

± 

0.000 

SN1991T 

21.454 

± 

0.008 

Pastorello et al. (2007) 

20.558 

± 

0.040 

SN2011fe 

19.001 

± 

0.027 

Wang et al. (2008b) 

19.196 

± 

0.040 


20.175 

± 

0.107 


20.410 

± 

0.123 
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Table B3 (cont'd) 


SN 

Phase t 
(days) 

Obs. Date 

m B (t) 

(mag) 

Phot. Ref. 

SN2006X 

360 

20070214 

20.824 ±0.151 


SN2007af 

103 

20070620 

17.181 ±0.005 

SN2011fe 

SN2007af 

108 

20070625 

17.247 ± 0.006 


SN2007af 

120 

20070707 

17.406 ± 0.007 


SN2007af 

123 

20070710 

17.446 ± 0.008 


SN2007af 

128 

20070715 

17.512 ±0.008 


SN2007af 

131 

20070718 

17.552 ± 0.009 


SN2007af 

151 

20070807 

17.816 ±0.012 


SN2007af 

165 

20070821 

18.002 ±0.014 


SN2007af 

308 

20080111 

19.894 ± 0.037 


SN2007gi 

161 

20080115 

17.920 ±0.014 

Zhang et al. (2010) 

SN20071e 

317 

20080827 

20.653 ± 0.000 

SN1991T 

SN2007sr 

177 

20080623 

17.994 ±0.021 

CfA, LOSS 

SN2009ig 

405 

20101017 

21.729 ±0.050 

CfA 

SN20091e 

324 

20101016 

22.651 ±0.050 

CfA 

SN2010ev 

178 

20101231 

20.199 ±0.050 

T15b 

SN2010ev 

272 

20110404 

21.183 ±0.050 


SN2010gp 

279 

20110501 

23.150 ±0.050 

T15b 

SN2010hg 

199 

20110403 

19.765 ± 0.050 

T15b 

SN2010kg 

289 

20110927 

22.592 ± 0.050 

T15b 

SN2011K 

341 

20111228 

23.096 ± 0.050 

T15b 

SN201 lae 

310 

20120101 

20.671 ±0.050 

T15b 

SN201 lat 

349 

20120227 

21.656 ±0.050 

T15b 

SN201 lby 

206 

20111202 

18.337 ± 0.022 

SN2011fe 

SN201 lby 

310 

20120315 

19.713 ±0.039 


SN201 lek 

423 

20121011 

24.328 ± 0.050 

T15b 

SN201 lfe 

074 

20111123 

13.630 ± 0.004 

Richmond & Smith (2012) 

SN201 lfe 

114 

20120102 

14.157 ±0.007 


SN201 lfe 

196 

20120324 

15.238 ±0.019 


SN201 lfe 

230 

20120427 

15.686 ± 0.024 


SN201 lfe 

276 

20120612 

16.292 ± 0.030 


SN201 lfe 

314 

20120720 

16.792 ± 0.036 


SN201 lim 

314 

20121016 

23.522 ± 0.050 

T15b 

SN201 liv 

318 

20121024 

20.717 ±0.050 

T15b 

SN201 ljh 

414 

20130221 

24.131 ±0.050 

T15b 

SN2012cg 

330 

20130507 

19.824 ± 0.085 

SN1994ae 

SN2012cg 

342 

20130513 

19.917 ±0.087 


SN2012cu 

340 

20130603 

23.060 ± 0.050 

T15b 

SN2012fr 

101 

20130221 

15.696 ±0.015 

Zhang et al. (2014) 

SN2012fr 

116 

20130308 

15.927 ±0.020 


SN2012& 

125 

20130317 

16.065 ± 0.024 


SN2012fr 

151 

20130412 

16.466 ± 0.034 


SN2012& 

222 

20130622 

17.560 ± 0.065 


SN2012fr 

261 

20130731 

18.161 ±0.082 


SN2012fr 

340 

20131018 

19.379 ±0.117 


SN2012fr 

357 

20131103 

19.409 ±0.100 


SN2012fr 

357 

20131103 

19.625 ±0.124 


SN2012fr 

367 

20131114 

19.795 ±0.129 


SN2012hr 

283 

20131006 

20.487 ± 0.005 

SN2003du 

SN2012hr 

368 

20131230 

21.749 ±0.007 


SN2012ht 

437 

20140313 

22.714 ±0.100 

T15b 

SN2013aa 

137 

20130710 

15.768 ± 0.000 

SN1991T 

SN2013aa 

185 

20130827 

16.438 ± 0.000 


SN2013aa 

202 

20130913 

16.676 ± 0.000 


SN2013aa 

342 

20140131 

18.633 ±0.001 


SN2013aa 

358 

20140216 

18.857 ±0.001 


SN2013aa 

430 

20140422 

19.765 ± 0.001 
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Table B3 (cont’d) 


SN 

Phase t 
(days) 

Obs. Date 

m B (t) 

(mag) 

Phot. Ref. 

SN2013cs 

320 

20140322 

20.909 ± 0.068 

SN2002fk 

SN2013cs 

322 

20140324 

20.939 ± 0.069 


SN2013cs 

322 

20140324 

21.475 ±0.100 


SN2013ct 

223 

20131218 

18.144 ±0.100 

T15b 

SN2013dl 

160 

20141130 

20.034 ±0.100 

T15b 

SN2013dl 

184 

20141224 

20.331 ±0.100 


SN2013dy 

333 

20140626 

20.716 ±0.033 

Pan et al. (2015a) 

SN2013dy 

419 

20140920 

22.064 ± 0.045 


SN2013ef 

174 

20141224 

21.511 ±0.100 

T15b 

SN2013gy 

276 

20140920 

21.242 ±0.062 

SN2007af 

SN2014J 

231 

20140920 

16.956 ±0.119 

Folev et al. (2014) 


Note. — For SNe with poor late-time photometric coverage, the late light curve 
of a surrogate SN (denoted in boldface in the Photometry Reference column) is em¬ 
ployed. _____ 


CfA: Riess et al. (19991); 

Jha et al. (2006); Hicken et al. (2009, 2012J); LOSS 

Ganeshalingam et al. (2010 

); CSP: Contreras et al. |2()ld); Stritzinger et al. (1201 lb 
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Table B4. Co line fluxes from spectra 


SN 

Date 



Co Line Flux 
(erg cm 2 s~ 1 ) 





B-Band Flux 
(erg cm 2 s —1 ) 



l^min 

(A) 

lOmax 

(A) 

SN1990N 

19901217 

6.06 

X 

IQ - 16 

±3.39 

X 

10-rs 

3.52 

X 

10 -15 

± 1.95 

X 

10- 17 

5710.5 

6038.4 

SN1990N 

19910112 

5.99 

X 

10- 14 

± 1.76 

X 

10- 17 

3.98 

X 

10- 13 

±7.37 

X 

10- 17 



SN1990N 

19910222 

3.69 

X 

10 -14 

±8.73 

X 

10“ 17 

3.30 

X 

10” 13 

± 5.29 

X 

IO -16 



SN1990N 

19910322 

9.78 

X 

10 -15 

± 4.25 

X 

10“ 17 

1.28 

X 

10- 13 

± 2.03 

X 

lO-re 



SN1990N 

19910416 

7.15 

X 

IO" 15 

±5.64 

X 

10- 17 

1.03 

X 

10- 13 

± 2.24 

X 

IO" 16 



SN1990N 

19910608 

1.89 

X 

10 -15 

±4.67 

X 

10- 17 

5.34 

X 

IO- 14 

± 1.46 

X 

IO -16 



SN1991M 

19910820 

1.67 

X 

10—i6 

±5.39 

X 

10- 18 

8.47 

X 

10—i6 

± 1.29 

X 

10~ 17 

5770.8 

6031.1 

SN1991T 

19910819 

4.99 

X 

10~ 14 

± 1.98 

X 

10-re 

1.28 

X 

10- 13 

± 5.48 

X 

lO-re 

5664.2 

6121.0 

SN1991T 

19911031 

1.34 

X 

10“ 15 

±2.23 

X 

10- 17 

8.20 

X 

10“ 15 

±4.80 

X 

10- 17 



SN1991T 

19920111 

5.49 

X 

10” 14 

±8.12 

X 

10“ 17 

4.57 

X 

10- 13 

± 2.40 

X 

10~ 16 



SN1991T 

19920313 

1.18 

X 

lO-m 

±4.46 

X 

10- 18 

1.09 

X 

IO" 15 

±6.92 

X 

IO” 18 



SN1991T 

19920411 

5.67 

X 

io - 17 

±2.51 

X 

10- 18 

8.32 

X 

10—i6 

±3.12 

X 

io - 17 



SN1993Z 

19940317 

2.17 

X 

io - 16 

±3.14 

X 

10- 18 

1.14 

X 

IO" 15 

± 1.00 

X 

io - 17 

5757.8 

6081.7 

SN1993Z 

19940418 

9.90 

X 

10- 17 

±2.47 

X 

10-rs 

9.73 

X 

10—i6 

±5.87 

X 

10-rs 



SN1994D 

19951124 

4.69 

X 

10” 15 

±3.19 

X 

10~ 16 

4.88 

X 

IQ - 14 

±4.08 

X 

10~ 16 

5732.7 

6070.7 

SN1994ae 

19950422 

6.74 

X 

10—i6 

±7.37 

X 

10- 18 

3.77 

X 

IO" 15 

±9.61 

X 

10- 18 

5747.5 

6027.4 

SN1994ae 

19950501 

5.38 

X 

IO" 14 

± 1.84 

X 

IO -16 

4.04 

X 

10- 13 

±6.40 

X 

IO" 16 



SN1995D 

19951124 

4.96 

X 

10” 15 

±3.49 

X 

10~ 16 

4.88 

X 

IO” 14 

±4.08 

X 

10~ 16 

5677.7 

6080.5 

SN1995D 

19951202 

1.96 

X 

10” 15 

± 1.70 

X 

10~ 16 

2.67 

X 

IQ - 14 

± 2.00 

X 

10~ 16 



SN1996X 

19970210 

2.45 

X 

IO" 15 

±3.74 

X 

10- 17 

2.97 

X 

IO" 14 

± 1.09 

X 

IO" 16 

5681.8 

6075.7 

SN1998aq 

19981124 

3.80 

X 

IO" 14 

± 1.20 

X 

IO' 15 

2.68 

X 

10- 13 

± 2.37 

X 

IO” 15 

5668.6 

6077.6 

SN1998aq 

19981214 

5.90 

X 

IO" 14 

±4.93 

X 

10~ 16 

3.66 

X 

10- 13 

± 1.01 

X 

10“ 15 



SN1998aq 

19981224 

4.86 

X 

IO” 14 

± 7.15 

X 

10~ 16 

3.86 

X 

10” 13 

± 1.79 

X 

10” 15 



SN1998bu 

19981114 

1.27 

X 

10- 13 

±7.07 

X 

IO -16 

5.36 

X 

10- 13 

± 1.95 

X 

IO” 15 

5706.7 

6055.0 

SN1998bu 

19981125 

1.39 

X 

10- 13 

±5.29 

X 

IO -16 

6.45 

X 

10- 13 

± 1.20 

X 

10- 15 



SN1998bu 

19981213 

8.76 

X 

IO” 14 

±3.94 

X 

10~ 16 

5.13 

X 

10- 13 

±9.64 

X 

10~ 16 



SN1998bu 

19981222 

8.04 

X 

IO” 14 

±6.98 

X 

10~ 16 

4.24 

X 

10- 13 

± 1.59 

X 

10“ 1S 



SN1998bu 

19990110 

4.58 

X 

10—i6 

± 1.19 

X 

10- 18 

3.71 

X 

IO" 15 

± 5.44 

X 

10- 18 



SN1998bu 

19990117 

5.47 

X 

IQ - 14 

±3.47 

X 

IO" 16 

2.95 

X 

10- 13 

±8.72 

X 

IO" 16 



SN1998bu 

19990223 

2.08 

X 

10—i6 

± 1.23 

X 

10- 18 

1.72 

X 

IO" 15 

±3.75 

X 

10- 18 



SN1998bu 

19990413 

6.52 

X 

10” 15 

± 7.18 

X 

10“ 17 

6.01 

X 

IQ - 14 

± 1.06 

X 

10~ 16 



SN1998bu 

19990424 

1.42 

X 

10- 16 

±2.54 

X 

10- 18 

1.04 

X 

IO" 15 

± 7.50 

X 

10- 18 



SN1999aa 

19991109 

1.80 

X 

10- 16 

± 1.12 

X 

10- 18 

8.69 

X 

10- 16 

± 2.46 

X 

10- 18 

5616.9 

6197.1 

SN1999aa 

19991205 

7.22 

X 

10- 17 

±8.10 

X 

IO" 19 

3.46 

X 

IO" 16 

± 2.46 

X 

10- 18 



SN2002cs 

20021106 

1.34 

X 

10- 16 

±6.21 

X 

10~ 19 

6.04 

X 

10- 16 

± 2.59 

X 

10- 17 

5655.5 

6097.2 

SN2002dj 

20030326 

3.41 

X 

10“ 15 

±2.40 

X 

10“ 17 

2.31 

X 

IO - 14 

±3.97 

X 

10~ 17 

5678.4 

6148.5 

SN2002er 

20030410 

5.02 

X 

IO" 15 

±8.81 

X 

10- 17 

4.18 

X 

IQ - 14 

± 2.10 

X 

IO" 16 

5709.1 

6086.8 

SN2002fk 

20030227 

4.95 

X 

10-i6 

±6.66 

X 

10- 19 

3.08 

X 

IO" 15 

± 2.47 

X 

10- 18 

5739.0 

6008.7 

SN2003du 

20030823 

1.36 

X 

10- 13 

±3.14 

X 

10~ 16 

5.08 

X 

10- 13 

± 1.06 

X 

10” 15 

5736.0 

6027.0 

SN2003du 

20030921 

4.90 

X 

IO - 14 

±9.36 

X 

10- 17 

2.47 

X 

10- 13 

±4.95 

X 

IO” 16 



SN2003du 

20030922 

4.24 

X 

IQ - 14 

± 1.16 

X 

IO" 16 

2.12 

X 

10- 13 

± 6.54 

X 

IO -16 



SN2003du 

20030925 

5.96 

X 

IO - 14 

± 1.36 

X 

lO-m 

3.00 

X 

10- 13 

± 7.96 

X 

IO' 16 



SN2003du 

20031201 

1.39 

X 

IO" 14 

± 1.19 

X 

lO-m 

1.05 

X 

10- 13 

±3.34 

X 

IO” 16 



SN2003du 

20031213 

8.91 

X 

IO” 15 

± 1.06 

X 

IO -16 

9.51 

X 

IO -14 

± 2.63 

X 

IO” 16 



SN2003du 

20040202 

3.34 

X 

IO" 15 

±6.14 

X 

10- 17 

6.68 

X 

IO" 14 

± 2.01 

X 

IO" 16 



SN2003du 

20040517 

5.05 

X 

10-i6 

± 7.93 

X 

10- 18 

1.24 

X 

IO -14 

± 1.96 

X 

10- 17 



SN2003gs 

20040214 

4.13 

X 

10- 17 

±2.53 

X 

10“ 19 

4.83 

X 

10- 16 

± 1.08 

X 

10“ 17 

5762.4 

6006.7 

SN2003hv 

20031228 

4.12 

X 

10- 13 

± 1.28 

X 

10“ 15 

1.64 

X 

10- 12 

± 2.68 

X 

10~ 15 

5687.6 

6061.5 

SN2003hv 

20040129 

3.33 

X 

10- 13 

± 1.31 

X 

IO” 15 

1.96 

X 

IO -12 

± 3.52 

X 

IO" 15 



SN2003hv 

20040725 

2.17 

X 

10~ 15 

± 1.56 

X 

10- 17 

3.00 

X 

IO - 14 

±3.21 

X 

10- 17 



SN2004bv 

20041114 

3.89 

X 

10- 16 

± 1.05 

X 

10- 18 

2.00 

X 

IO" 15 

± 1.87 

X 

10- 18 

5744.3 

6057.6 

SN2004eo 

20050516 

1.46 

X 

10 -15 

±4.14 

X 

10- 18 

1.22 

X 

IO - 14 

± 1.10 

X 

10- 17 

5737.5 

6034.0 

SN2005cf 

20060427 

4.44 

X 

10-rs 

±2.06 

X 

10- 19 

5.34 

X 

10- 17 

± 1.99 

X 

10- 19 

5764.7 

5992.4 

SN2006X 

20060626 

7.25 

X 

10-i6 

± 1.12 

X 

10- 17 

8.76 

X 

10-i6 

± 1.13 

X 

10- 17 

5670.0 

6118.6 

SN2006X 

20060721 

3.12 

X 

IO” 14 

±4.78 

X 

10~ 16 

5.35 

X 

IO” 14 

± 1.20 

X 

10“ 15 



SN2006X 

20061123 

7.25 

X 

10- 17 

±5.51 

X 

10“ 19 

2.01 

X 

10-i6 

± 7.62 

X 

10“ 19 



SN2006X 

20061223 

3.70 

X 

10- 16 

±2.24 

X 

10- 18 

5.19 

X 

10- 16 

± 5.02 

X 

10- 17 
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Table B4 (cont'd) 


SN 

Date 



Co Line Flux 
(erg cm 2 s —1 ) 





B-Band Flux 
(erg cm 2 s — *) 



'Wmin 

(A) 

'LL) max 

(A) 

SN2006X 

20070214 

5.51 

X 

10- 17 

± 1.42 

X 

10- 19 

1.50 

X 

10- 16 

±4.34 

X 

10- 19 



SN2007af 

20070620 

2.96 

X 

10“ 13 

± 6.13 

X 

10— 16 

1.09 

X 

IO -12 

± 1.63 

X 

IO" 15 

5762.2 

6027.6 

SN2007af 

20070625 

2.15 

X 

10 -13 

± 5.45 

X 

10- 16 

7.80 

X 

10- 13 

±2.19 

X 

IO" 15 



SN2007af 

20070707 

3.06 

X 

10- 15 

± 4.33 

X 

10 —i8 

1.36 

X 

IO" 14 

± 1.73 

X 

10- 17 



SN2007af 

20070710 

8.47 

X 

10" 14 

± 4.97 

X 

10- 16 

3.08 

X 

io- 13 

± 1.39 

X 

IQ - 15 



SN2007af 

20070715 

8.25 

X 

10 —i6 

± 2.73 

X 

10- 18 

4.00 

X 

IO" 15 

±9.48 

X 

10- 18 



SN2007af 

20070718 

7.49 

X 

10- 14 

± 6.87 

X 

10— 16 

3.47 

X 

10- 13 

±3.23 

X 

10~ 15 



SN2007af 

20070807 

2.50 

X 

IO” 16 

± 3.72 

X 

10 —i8 

1.67 

X 

10- 15 

± 1.04 

X 

10- 17 



SN2007af 

20070821 

3.57 

X 

IO" 16 

± 4.16 

X 

10- 18 

2.15 

X 

IO" 15 

± 1.93 

X 

10- 17 



SN2007af 

20080111 

2.19 

X 

10 -15 

± 4.58 

X 

10- 17 

3.08 

X 

IO" 14 

±8.19 

X 

10- 17 



SN2007gi 

20080115 

1.20 

X 

10~ 15 

± 3.12 

X 

10- 18 

5.08 

X 

IO" 15 

±8.44 

X 

IO” 18 

5706.8 

6146.6 

SN20071e 

20080827 

1.80 

X 

10- 17 

± 4.41 

X 

io~ 19 

2.02 

X 

10—i6 

±5.68 

X 

IO- 19 

5721.8 

6087.5 

SN2007sr 

20080623 

4.42 

X 

10- 14 

± 5.84 

X 

10" 17 

3.05 

X 

IO -13 

±4.04 

X 

IO" 16 

5765.7 

6059.7 

SN2008Q 

20080828 

9.20 

X 

10- 17 

± 2.08 

X 

io- 19 

6.25 

X 

10—i6 

±5.91 

X 

io- 19 

5691.4 

6089.5 

SN2009ig 

20101017 

1.50 

X 

10- 15 

± 6.87 

X 

10- 18 

1.23 

X 

IO" 14 

± 1.85 

X 

io- 17 

5496.8 

6156.7 

SN20091e 

20101016 

5.18 

X 

IO” 16 

± 9.02 

X 

10- 18 

5.27 

X 

IO" 15 

±2.33 

X 

io- 17 

5587.7 

6195.9 

SN2010ev 

20101231 

1.23 

X 

10- 14 

± 2.27 

X 

io- 17 

5.05 

X 

IO" 14 

±5.09 

X 

io- 17 

5649.8 

6152.9 

SN2010ev 

20110404 

3.21 

X 

10 -15 

± 8.29 

X 

10- 18 

2.04 

X 

IO" 14 

±3.81 

X 

io- 17 



SN2010gp 

20110501 

2.82 

X 

IO” 16 

± 4.23 

X 

10- 18 

3.33 

X 

IO” 15 

± 1.18 

X 

10“ 17 

5726.6 

6063.4 

SN2010hg 

20110403 

9.34 

X 

10- 15 

± 2.08 

X 

io- 17 

7.52 

X 

IO" 14 

± 7.41 

X 

io- 17 

5728.8 

6037.9 

SN2010kg 

20110927 

7.55 

X 

10 —i6 

± 9.49 

X 

10- 18 

5.57 

X 

IO" 15 

±3.43 

X 

io- 17 

5650.3 

6165.6 

SN201 lae 

20120101 

2.70 

X 

10 -15 

± 2.21 

X 

io- 17 

3.27 

X 

IO" 14 

±7.28 

X 

io- 17 

5684.3 

6087.6 

SN2011at 

20120227 

1.93 

X 

IO" 15 

± 8.90 

X 

10- 18 

1.32 

X 

IO- 14 

±2.58 

X 

io- 17 

5618.1 

6118.7 

SN2011by 

20111202 

1.85 

X 

10 -15 

± 4.98 

X 

10- 18 

1.53 

X 

IO" 14 

± 1.68 

X 

10“ 17 

5723.6 

6039.5 

SN201 lby 

20120315 

5.84 

X 

10- 17 

± 4.32 

X 

io- 19 

8.72 

X 

IO" 16 

± 1.26 

X 

10- 18 



SN2011ek 

20121011 

1.73 

X 

10 —i6 

± 4.75 

X 

10- 18 

1.12 

X 

IO" 15 

±9.34 

X 

10- 18 

5544.0 

6130.4 

SN2011fe 

20111123 

7.06 

X 

10- 12 

± 1.78 

X 

10 _ls 

1.76 

X 

IO- 11 

±5.31 

X 

10 _ls 

5738.3 

6027.2 

SN2011fe 

20120102 

1.38 

X 

10- 12 

± 2.86 

X 

IO -16 

5.87 

X 

10~ 12 

± 1.76 

X 

10~ 15 



SN201 lfe 

20120324 

2.14 

X 

10- 13 

± 5.88 

X 

10- 17 

1.81 

X 

10- 12 

±5.28 

X 

IO” 16 



SN201 lfe 

20120427 

8.30 

X 

10- 14 

± 2.34 

X 

10- 17 

8.94 

X 

10- 13 

±2.50 

X 

IO" 16 



SN2011fe 

20120612 

5.69 

X 

10- 14 

± 1.77 

X 

10- 17 

8.42 

X 

10- 13 

±2.62 

X 

10— 18 



SN2011fe 

20120720 

2.83 

X 

10- 14 

± 4.42 

X 

10“ 17 

5.26 

X 

10- 13 

±3.84 

X 

10” 16 



SN201 lim 

20121016 

2.80 

X 

IO" 16 

± 7.58 

X 

10- 18 

2.36 

X 

10 -15 

± 1.93 

X 

10- 17 

5687.0 

6119.4 

SN2011iv 

20121024 

2.77 

X 

10 -15 

± 4.00 

X 

10- 17 

3.13 

X 

IO" 14 

±8.86 

X 

10- 17 

5693.8 

6089.2 

SN2012cg 

20130507 

3.98 

X 

10 -15 

± 5.08 

X 

10- 17 

5.30 

X 

IO- 14 

±6.44 

X 

10- 17 

5698.9 

6113.3 

SN2012cg 

20130513 

4.30 

X 

10- 15 

± 3.59 

X 

10- 17 

6.51 

X 

IO" 14 

± 1.22 

X 

IO -16 



SN2012cu 

20130603 

6.26 

X 

10-i6 

± 7.42 

X 

10- 18 

3.62 

X 

IO" 15 

±2.06 

X 

IO” 17 

5692.6 

6096.8 

SN2012fr 

20130221 

6.02 

X 

IO -13 

± 6.43 

X 

IO -16 

1.97 

X 

IO -12 

±2.79 

X 

IO" 15 

5759.2 

6061.5 

SN2012fr 

20130308 

3.30 

X 

10- 13 

± 4.81 

X 

IO -16 

1.28 

X 

10- 12 

±2.34 

X 

10~ ls 



SN2012fr 

20130317 

3.52 

X 

10“ 13 

± 4.66 

X 

IO” 16 

1.42 

X 

10- 12 

± 1.48 

X 

IO” 15 



SN2012fr 

20130412 

2.22 

X 

10- 13 

± 1.51 

X 

10 -15 

1.49 

X 

10- 12 

±4.85 

X 

10 -15 



SN2012fr 

20130622 

3.52 

X 

10- 14 

± 3.92 

X 

10— 16 

3.57 

X 

IO -13 

± 1.62 

X 

10~ 15 



SN2012fr 

20130731 

1.31 

X 

10- 14 

± 1.56 

X 

10- 16 

1.77 

X 

10- 13 

±5.51 

X 

10- 16 



SN2012fr 

20131018 

3.28 

X 

10- 15 

± 1.08 

X 

IO -16 

6.72 

X 

IO" 14 

±5.55 

X 

IO" 16 



SN2012fr 

20131103 

5.47 

X 

10~ 15 

± 2.94 

X 

10- 18 

1.19 

X 

10- 13 

± 1.98 

X 

10- 17 



SN2012fr 

20131114 

1.25 

X 

10 -15 

± 1.05 

X 

IO" 16 

4.16 

X 

IO" 14 

±3.91 

X 

IO" 16 



SN2012hr 

20131006 

3.61 

X 

10- 17 

± 2.28 

X 

10- 19 

5.34 

X 

IO" 16 

± 7.39 

X 

10- 19 

5775.3 

6041.1 

SN2012hr 

20131230 

2.53 

X 

IO” 16 

± 2.76 

X 

10” 17 

7.26 

X 

IO" 15 

± 1.39 

X 

10” 16 



SN2013aa 

20130710 

3.75 

X 

10- 13 

± 4.69 

X 

IO -16 

1.60 

X 

10~ 12 

± 1.07 

X 

10~ 15 

5720.3 

6055.9 

SN2013aa 

20130827 

5.88 

X 

10- 14 

± 1.19 

X 

10~ 16 

4.14 

X 

10- 13 

±5.40 

X 

IO" 16 



SN2013aa 

20130913 

6.54 

X 

10- 14 

± 4.35 

X 

IO" 16 

5.13 

X 

10- 13 

± 1.62 

X 

IO" 15 



SN2013aa 

20140131 

1.76 

X 

10- 15 

± 1.48 

X 

10— 16 

5.26 

X 

IO" 14 

±3.92 

X 

IO" 16 



SN2013aa 

20140216 

2.10 

X 

10~ 15 

± 2.73 

X 

10- 17 

3.78 

X 

IO" 14 

±4.23 

X 

10” 17 



SN2013aa 

20140422 

8.32 

X 

10 -16 

± 2.73 

X 

10- 17 

1.61 

X 

IO” 14 

±3.68 

X 

10- 17 



SN2013cs 

20140322 

3.28 

X 

10 —i6 

± 9.26 

X 

IO" 17 

7.49 

X 

IO" 15 

±3.94 

X 

IO" 16 

5735.8 

6061.7 

SN2013cs 

20140324 

9.99 

X 

IO” 16 

± 1.80 

X 

10- 17 

1.56 

X 

IO -14 

±2.24 

X 

10“ 17 



SN2013ct 

20131218 

3.71 

X 

10- 14 

± 8.54 

X 

10- 17 

3.35 

X 

10- 13 

±2.32 

X 

IO" 16 

5753.1 

6030.3 
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Table B4 (cont'd) 


SN 

Date 

Co Line Flux 
(erg cm 2 s —1 ) 

B-Band Flux 
(erg cm 2 s —1 ) 

(A) 

1C! max 

(A) 

SN2013dl 

20141130 

1.17 x 10~ 14 ±6.49 x 10- 17 

5.87 x 10 -14 ± 1.04 x 10~ 16 

5763.5 

6061.3 

SN2013dl 

20141224 

7.29 x 10“ 15 ± 6.95 x 10 -17 

4.47 x 10” 14 ± 8.46 x 10 -17 



SN2013dy 

20140626 

5.54 x lO" 15 ± 3.61 x 10~ 17 

6.29 x 10 -14 ± 3.65 x 10“ 15 

5636.9 

6168.7 

SN2013dy 

20140920 

1.45 x 10“ 14 ± 1.34 x 10~ 16 

1.00 x 10- 13 ± 6.47 x 10- 15 



SN2013ef 

20141224 

3.28 x 10 -15 ± 3.58 x 10 -17 

1.51 x 10 -14 ± 5.48 x 10~ 17 

5733.8 

6053.0 

SN2013gy 

20140920 

7.63 x 10 -15 ± 5.97 x 10~ 17 

5.72 x 10 -14 ± 1.82 x 10- 15 

5726.8 

6064.7 

SN2014J 

20140920 

1.65 x 10” 12 ± 7.92 x 10- 15 

5.19 x 10 -12 ±9.56 x 10- 15 

5737.8 

6069.5 
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